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INTRODUCTION 


'  y  There  Is  currently  considerable  interest  In  the  use  of  nitramine 
oxidizers  in  gun  and  rocket  propellants 'Tsee  for  example  Reference  1) ,  'due  to 
their  desirable  properties,  such  as  absence  of  noxious  combustion  products, 
high  specific  impulse  and  impetus^,  and  thermal  stability.  The  most  commonly 
used  nitramine  oxidizers  are  RDXl I) , (also  known  by  such  names  as  hexogen  or 
cyclonite) and -currently  indexed  in  Chemical  Abstracts  under  the  heading 
"l,3,5-Triazini7~~hexahydro-l,  3^  5-trinitro"  and^HMX  (II),  also  known  by  such 
names  as  Octogen  and  indexed  in  Chemical  Abstracts  finder  the  heading  ”1,3, 5, 7- 
.Tetrazocine,  octahydro-1 ,3,5,7-tetranitro."  -ifoc  7*  ^ )  --•>  .  H  ^ 

Before  the  combustion  snd  deflagration  behavior  of  these  materials  can  be 
understood  and  intelligently  modified,  it  will  be  necessary  to  have  an 
understanding  of  their  thermal  decomposition  chemistry.  This  is  true  for 
several  reasons.  First,  combustion  modellers  use  the  kinetics  and  kinetic 
parameters  as  input  into  their  models.  Second,  the  product  distributions  from 
the  decomposition  reactions  are  also  needed  as  input  for  many  of  these 
models.  Also,  an  understanding  of  the  detailed  chemical  decomposition  modes 
might  well  make  it  possible  to  suggest  new  types  of  additives  for  combustion 
modification. 

At  present,  the  chemical  mechanisms  involved  in  the  decomposition  of 
these  materials  are  not  at  all  well  understood.  There  are  a  number  of 
apparent  conflicts  in  the  literature  with  regard  to  such  topics  as  kinetic 
parameters,  product  distributions,  autocatalysis  and  acceleration  and/or 
inhibition  of  the  decomposition  by  substances,  some  cf  which  are  known 
products  of  the  decomposition.  Many  of  these  apparent  conflicts  probably 
arise  from  the  great  complexity  of  the  decomposition.  The  course  of  the 
reaction  seems  to  be  very  dependent  on  experimental  conditions  due  to  the 
apparent  operation  of  simultaneous  gas,  liquid  and  solid-phase  decomposition, 
leading  to  products  many  of  which  are  able  to  react  among  themselves  and  also 
to  accelerate  or  inhibit  the  decomposition  of  starting  material.  In  addition, 
the  nature  of  even  the  first  step  of  the  decomposition  reaction  is  not  at  all 
well  understood,  quite  possibly  as  a  result  of  variations  with  temperature, 
pressure  and  state  of  aggregation  (gas,  liquid,  solid,  dissolved  in  solvent) 
in  the  relative  rates  of  the  many  reactions  (N-NO2  cleavage,  HONO  elimination, 
various  forms  of  C-N  cleavage,  etc.)  that  have  been  proposed. 

In  view  of  these  considerations,  it  appeared  that  there  was  a  need  for  a 
critical,  interpretative  review  of  the  literature  on  decomposition  of 
nltramlnes,  with  the  emphasis  on  HMX  and  RDX.  Such  a  review  would  summarize 
what  is  known  about  the  chemical  and  physical  mechanisms  involved,  and 
involves  an  attempt  to  resolve  conflicts  and  to  suggest  new  experiments  to 
elucidate  the  mechanisms  involved. 

The  present  report  is  one  of  a  series  describing  such  a  review.  This 
review  is  written  from  the  point  of  view  of  an  organic  chemist  and  so  is 
complementary  to  reviews  such  as  that  by  McCarty,  which  emphasize  combustion 


1K.  P.  McCarty,  "HMX  Propellant  Combustion  Studies,"  AFRPL-TR-76-59  (AD-BO  17 
527L). 
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behavior  as  seen  by  a  physical  chemist  or  combustion  engineer.  In  the  present 
report,  available  information  on  kinetic  parameters  for  HMX  and  RDX 
decomposition  is  summarized,  and  an  attempt  is  made  to  arrive  at  suggested 
values  for  activation  energies  and  frequency  factors,  and  to  discuss  the 
results  in  terms  of  the  possible  chemical  and  physical  mechanisms  involved  in 
the  decomposition  reactions.  Brief  preliminary  discussions  of  possible 
chemical  mechanisms  and  of  needed  research  have  already  appeared, 2  and 
subsequent  reports  in  the  present  series  will  deal  with  such  topics  as  product 
distributions  and  autocatalysis.  There  have  been  a  number  of  reviews  of  HMX 
and  RDX  decomposition,-^-^  some  of  which  came  to  our  attention  after  the 
present  work  was  in  progress,  but  there  still  seems  to  be  a  definite  need  for 
a  critical,  interpretative  review  written  from  the  point  of  view  of  an  organic 
or  physical-organic  chemist. 

In  the  present  report,  an  attempt  has  been  made  to  provide  comprehensive 
literature  coverage  to  approximately  mid-1980,  although  some  later  reports  and 
papers  are  included. 

The  scope  and  organization  of  the  present  report  are  as  follows:  The 
available  values  for  activation  energy  and  frequency  factor  A  (expressed  as 
log  A)  are  summarized  in  Appendices  A-D.  Appendix  A  summarizes  the  very  few 
available  values  for  vapor-phase  decomposition  of  HMX  and  RDX.  Append!::  B 
summarizes  data  for  HMX  decomposition  in  the  neat  liquid  and  dissolved 
states.  Appendix  C  similarly  summarizes  data  for  RDX  decomposition  in  the 
neat  liquid  and  dissolved  states.  Appendix  D  covers  decomposition  of  HMX  and 
RDX  at  temperatures  below  their  respective  melting  points.  These  Appendices 


^(a)  M.  A.  Schroeder,  "Critical  Analysis  of  Nitramine  Decomposition  Results: 
Some  Comments  on  Chemical  Mechanisms,"  Proceedings,  16th  JANNAF  Combustion 
Meeting,  Naval  Postgraduate  School,  Monterey,  CA,  Sep  10-14,  1979,  Voi.  2, 
p.  17;  (b)  M.  A.  Schroeder,  "Critical  Analysis  of  Nitramine  Decomposition 
Data:  Some  Suggestions  for  Needed  Research  Work,”  Memorandum  Report  ARBRL- 
MR-  3181,  Jure  1982  (AD-A116  194). 

-*J.  Wilby,  "Thermal  Decomposition  of  RDX/TNT  Mixtures,  Part  1,"  A.R.D.E. 
Report  (MX)  16/59,  August  1959  (AD-285  991). 

*F.  C.  Rauch  and  R.  B.  Wainright,  "Studies  on  Composition  B,"  Final  Report 
Contract  No.  DAAA21-68-C-0334,  American  Cyanamid  Com' any,  Feb.  1969  (AD-850 
928). 

-*M.  Benreuven  and  L.  H.  Caveny,  "Nitramine  Monopropellant  Deflagration  and 
General  Nonsteady  Reacting  Rocket  Chamber  Flows,"  MAE  Report  No.  1455, 
Princeton  University,  Princeton,  NJ,  January  1980. 

**Ye.  Yu  Orlova,  N.  A.  Orlova,  et  al.,  "Octogen-Thermoresistant  Explosive," 
Publishing  House  "Nedra,"  Moscow,  1975;  FTD-ID  (RS)  T-0667-80,  2  May  1980 
(AD-B047  181). 
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are  seoarately  accompanied  by  discussions  Intended  to  justify  the  selection  of 
best  values  and  weights  of  data  from  each  title.  The  ’ference  numbers  in 
these  discussions  correspond  to  references  given  In  thi.c  particular 
Appendix.  A  more  general  discussion  of  sources  of  error,  chemical  mechanism, 
etc.  is  given  in  the  main  report;  this  more  general  discussion  Is  Intended  to 
tie  together  the  Information  in  ail  of  the  tables,  and  to  present  suggested 
values  for  kinetic  parameters  of  hMX  and  RDX  decomposition  in  the  liquid  and 
gaseous  states. 

It  should  be  understood  that  the  entire  decomposition  process  of  HMX 
cannot  be  described  bv  a  single  set  of  kinetics  constants;  the  kinetics  and 
parameters  differ  appreciably  with  experimental  conditions  such  as 
confinement,  temperature,  pressure  and  phase  (solid,  liquid,  gas,  solution) 
and  with  extent  of  reaction.  This  is  exFcclally  true  for  the  physically- 
complex  below-the-melting-point  decomposition.  In  which  considerable  evidence 
seems  consistent  with  the  idea  that  initially,  vaporization  and  gas-phase 
decomposition  (In  at  least  some  cases  probably  proceeding  side-by-side  with 
true  solid-state  decomposition)  leads  to  formation  of  liquid  products  which 
condense  onto  undecomposed  solid  and  induce  liquefaction,  leading  to 
acceleration  of  condensed-phase  decomposition.  However  the  rates  differ 
appreciably  even  between  the  homogeneous  solid,  liquid,  vapor  and  solution 
phases,  and  in  many  cases  even  the  rate  constants  for  homogeneous,  single- 
phase  decomposition  seems  subject  to  some  autocatalysis.  Thus,  wnen  rate  data 
such  as  those  summarized  here  are  used  for  predictive  purposes,  care  should  be 
taken  to  understand  as  well  as  possible  the  physical  and  chemical  mechanisms 
of  HMX  and  RDX  decomposition  under  the  conditions  of  pressure  and  temperature 
for  which  predictions  are  being  made;  and  to  take  as  much  account  as  possible 
of  variations  in  physical  and  chemical  mechanisms  with  temperature  and 
pressure. 

| 

It  should  be  remembered  that  the  activation  energies  and  frequency 
factors  given  here  are  the  result  of  efforts  to  arrive  at  values 
characteristic  of  certain  specific  chemical  decomposition  reactions  occurring 
in  the  homogeneous  liquid  or  vapor  {phases  (but  see  the  discussion  in  the 
preceding  paragraph).  Consequently,  they  are  strictly  applicable  to  complex 
processes  (such  as  combustion)  only  to  the  extent  that  these  processes  are 
representable  in  terms  of  these  same  specific,  homogeneous  decomposition 
steps.  Therefore,  when  the  values  jgiven  here  are  used  as  input  for  combustion 
or  explosive  modeling  studies,  they  can  be  expected  to  prove  most  useful  when 
used  in  conjunction  with  other  measured  or  estimated  parameters,  in  models 
which  take  detailed  account  of  the  individual  chemical  and  physical  processes 
involved.  There  is  a  need  for  further  development  of  propellant-combustion 
and  explosive  models  which  take  detailed  account  of  the  chemical  decomposition 
processes  involved,  especially  since  relative  rates  of  various  initial  and 
follow-up  steps  may  change  with  changes  in  pressure  and  temperature. 

II.  SOURCES  OF  ERROR 

One  important  source  of  uncertainty  in  studies  such  as  those  summarized 
here  is  that  use  of  a  very  short  temperature  range  for  the  Arrhenius  plot  may 
lead  to  inexact  values  for  the  frequency  factor  and  activation  energies. 

Other  important  sources  of  error  Include  self-heating;  vaporization  followed 
by  vapor  phase  decomposition  or  by  condensation  onto  the  cooler  parts  of  the 
apparatus;  reactions  taking  place  on  or  catalyzed  by  the  container  surface; 
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and  autocatalysis  or  autoinhibition.  Another  problem  is  that  it  is  ofter  not 
stated  whether  or  what  sort  of  statistical  analysis  was  employed.  This  can 
make  a  difference,  as  is  shown  by  the  following  example:  an  activation  energy 
of  48.7  kcal/mole  is  reported^ for  the  thermal  decomposition  of  neat  liquid 
RDX.  In  a  later  report”  in  the  same  series,  it  was  stated  thaw  least  squares 
analysis  of  the  data  gave  a  log  A  (sec~*)  of  22.4  and  an  activation  energy  of 
55.6  kcal/mole.  This  is  a  significant  change  from  the  earlier  values  which 
were  stated  to  be  the  result  of  a  preliminary  analysis.  This  example 
illustrates  the  possible  effect  of  statistical  treatments  on  the  results. 


A.  Temperature  Range 

Activation  energy  and  frequency  factors  are  determined  from  a  logari •;  ..tic 
plot  of  rate  constant  k  vs  reciprocal  temperature;  clearly  the  longer  the 
temperature  range  the  more  reliable  will  be  the  extrapolation  t'*  the  much 
higher  temperatures  characteristic  of  combustion  and  explosive  behavior.  Most 
authors  studying  HMX  and  RDX  decomposition  kinetics  have  worked  over 
'■elatively  short  temperature  ranges,  often  20°  or  less,  possibly  because 
temperature  dependence  due  to  the  relatively  high  (ca  50  kcal/mole)  activation 
energies  results  in  shortening  of  the  temperature  range  over  which  the 
decomposition  proceeds  at  conveniently  measurable  rates.  This  problem  is 
especially  severe  for  liquid  HMX,  which  is  already  within  its  decomposition 
range  at  the  melting  point  of  solid  HMX  (ca  280°).  In  spite  of  this  problem, 
Robertson^  was  able  to  study  neat  liquid  RDX  over  the  range  of  213-99°C  (AT  * 
86°)  and  HMX  over  the  range  271-314°C  ( AT  ■*  43°);  the  reaction  .t  higher 
temperatures  was  followed  using  a  membrane  manometer  attached  to  an  optical 
lever,  pressure-time  traces  being  recorded  on  a  moving  strip  of  photographic 
paper. 

In  the  present  report,  we  will  plot  the  most  reliable  data  on  common 
axes,  hoping  thereby  to  moot  the  question  of  differences  in  reliability  due  to 
use  of  different  temperature  ranges  in  the  different  studies  by  doing  common 
statistical  analysis  on  all  data  for  the  same  compound  and  phase. 


F.  C.  Rauch  and  A.  J.  Fanelli,  "The  Thermal  Decomposition  Kinetics  of 
Hexahydro-1 ,3 ,5-*‘rinitro-8-triazine  above  the  Melting  Point:  Evidence  for 
both  a  Gas  and  Liquid  Phase  Decomposition,”  J.  Phys.  Chem.,  Vol.  73, 
p.  1604,  1969. 

O 

F.  C.  Rauch  and  W.  P.  Colraan,  "Studies  on  Composition  B,”  Final  Report, 
Contract  No.  DAAA2 1-68-G-0334 ,  American  Cyanamid  Company,  March  1970  (AI>-869- 
226). 

a 

A.  J.  B.  Robertson,  "The  Thermal  Decomposition  of  Explosives,  Part  II. 
Cyclotrimethylenetrinitramine  and  Cyclotetraraethylenetetramitramine ,"  Trans . 
Faraday  Soc.,  Vol.  45,  p.  85,  1949. 
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B.  Self-Heating 


Self-heating  results  when  an  exothermic  reaction  generates  heat  faster 
than  the  surroundings  can  conduct  It  away.  Self-heating  would  be  favored^  by 
large  sample  sizes  and  by  higher  pressure  of  gas  over  the  sample  since  high 
pressure  would  lead  to  smaller  bubbles  of  decomposition  gases,  hence  to  less 
efficient  stirring  of  the  sample.  Thus  two  checks  for  self-heating  would  he 
whether  or  not  the  rate  were  Independent  of  (a)  sample  size  and  of  (b) 
pressure  of  Inert  gas  atmosphere.  Robertson9  found  no  self-heating  by  either 
methods  for  RDX  with  sample  sizes  from  4-45  mg,  or  for  HMX  with  sample  sizes 
from  1-7  mg.  Thus,  it  would  appear  that  in  manoraetric  kinetic  apparatus  of 
this  general  type,  use  of  samples  of  these  sizes  or  smaller  should  obviate  the 
problem  of  self-heating.  However,  as  Indicated  in  Appendices  B  and  C,  use  of 
larger  samples  did  sometimes  give  evidence  of  self-heating.  Furthermore, 
Goshgarlan  ^  in  describing  a  DSC  experiment  states  that  "sample  weight  (of 
HMX)  greater  than  1 .0  mg  increased  the  isothermal  temperature  during  HMX 
decomposition,  apparently  due  to  excess  thermal  evolution  beyond  the 
capabilities  of  the  instrument  to  adequately  regulate  a  constant  temperature." 
Thus  the  permissible  sample  size  in  a  DSC  apparatus  appears  to  be  much  smaller 
than  in  a  standard  kinetic  apparatus. 

If  self-heating  were  important,  it  should  give  artificially  high  rate 
constants  at  high  temperatures  resulting  In  high  values  for  activation  energy 
and  frequency  factor. 

C.  Sublimation  and  Ga3-Phase  Reactions 

Another  possible  important  source  of  error  in  condensed-phase  studies  Is 
vaporization  of  the  sample  followed  by  either  gas-phase  decomposition,  or  even 
sublimation  of  undecomposed  HMX  or  RDX  onto  cooler  portions  of  the  apparatus. 

In  most  of  the  kinetic  studies,  sublimation  was  retarded  either  by  carrying 
out  the  reaction  under  a  pressure  of  Inert  gas  or  by  carrying  out  the  reaction 
In  an  apparatus  designed  In  such  a  way  that  all  parts  were  evenly  heated  to 
the  temperature  of  the  experiments.  Thes**  precautions  serve  a  purpose,  since 
It  Is  widely  reported  that  in  vacuum  both  HMX  and  RDX  readily  vaporize  and 
condense  on  the  cooler  parts  of  the  apparatus  without  appreciable 
decomposition;  see  for  example  Robertson.’ 

The  vapor-phase  reaction  can  interfere  appreciably  with  determination  of 
liquid-phase  kinetic  parameters,  as  shown  by  the  results  of  Rogers,^  who 
found  that  the  vapor  phase  caused  a  quite  noticeable  hump  on  the  tall  of  the 


*^A.  J.  B.  Robertson,  "The  Thermal  Decomposition  of  Explosives,  Part  I. 
Ethylenedl nltramine  and  Tetryl,  "Trans.  Faraday  Soc.,  Vol.  44.  p.  677, 
1948.  '  “  ‘ 

B.  B.  Goshgarlan,  "The  Thermal  Decomposition  of  Cyclotrlmethylene- 
trlnitramine  (RDX)  and  Cyclotetramethylenetetranltramlne  (HMX),"  AFRPL-TR- 
78-76  ( AD- BO 32-27 5L). 

i  2 

R.  N.  Rogers,  "The  Determination  of  Condensed  Phase  Kinetics  Constants," 
Thermochtmlca  Acta,  Vol.  9,  p.444,  1974. 
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DSC  curve  of  decomposing  RDX.  The  gas-phase  material  disappears  Immediately 
after  the  liquid-phase  is  all  gone;  the  result  Is  a  hump  on  the  tail  of  the 
DSC  curve  at  the  point  where  liquid-phase  decomposition  ends  and  the  gas-phase 
RDX  starts  to  disappear.  While  liquid-phase  is  present,  what  is  actually 
being  observed  is  presumably  a  superposition  of  a  first-order  liquid  phase 
reaction  and  a  gas-phase  reaction  which  reduces  to  apparent  zero-order  due  to 
constant  volume  and  to  replenishment  of  gaseous  material  by  vaporization  of 
liquid.  In  the  case  of  RDX  this  led*^  to  an  apparent  Ea  »  43.1  kcal/mole  and 
A  =  2.44  x  10*®  sec-*.  When  the  vapor-phase  reaction  was  taken  account  of  by 
using  a  baseline  based  on  the  constant  gas-phase  decomposition,  Rogers 
obtained  Ea  =  47.1  kcal/mole  and  A  =*  2.02  x  10*®  sec-*,  in  good  agreement  with 
Robertson's  manometric  values.  In  spite  of  Roger's  publication,  subsequent 
workers  have  generally  continued  to  ignore  the  vapor-phase  contribution  to  the 
decomposition  of  neat  liquid  HMX  and  RDX. 

The  above  discussion  refers  to  isothermal  DSC.  When  temperature- 
programmed  nonisothermal  DSC  is  employed,  the  problem  becomes  even  worse  due 
to  probable  temperature  dependence  of  the  relative  importance  of  va;  *:-and 
liquid-phase  decomposition,  as  well  as  error  introduced  by  the  particular 
method  of  calculation  employed. 


Thus,  out  of  the  multitude  of  thermal  analysis  studies  of  Arrhenius 
parameters  of  HMX/RDX  decomposition,  the  gas-decoraposition-corrected  values  of 
Rogers*^***®  are  preferred.  The  4  kcal/mole  difference  between  Rogers' 
corrected  and  uncorrected  activation  energies  for  RDX  suggests  that  a  similar 
correction  should  be  applied  to  all  Isothermal  DSC  results;  examination  of 
Appendix  C  suggests  that  application  of  such  a  correction  would  bring  these 
values  into  the  same  range  as  those  found  from  non-DSC  kinetic  studies  on 
RDX.  See  also  the  discussion  in  Section  II.  F. 


The  effect  of  vapor-phase  decomposition  may  be  smaller  for  non-DSC 
studies.  Rauch  and  Wainwright  (see  page  29  of  Reference  4),  in  the  course  of 
their  reactant-disappearance  kinetic  studies  of  the  decomposition  of  liquid 
RDX,  calculated  (from  vapor  pressure)  the  percent  RDX  initially  in  the  gas 
phase  for  several  different  combinations  of  reactor  volume  and  initial  weight 
of  RDX  at  212°C.  These  percentages  varied  from  0.3  to  3.6,  and  seemed 
unrelated  to  the  rate  constant  which  varied  in  the  range  2.0±0.3  x  10-®  sec-* 
in  a  manner  unrelated  to  the  amount  of  RDX  present  in  the  gas  phase.  This  is 
surprising,  since  the  li  lid-phase  rate  constant  at  this  temperature  is 
calculated  to  be  1.25  x  lu-®  sec-*,  using  Robertson's  values’  (A  *  3.16  x 
10*®,  Ea  =»  47.5  kcal/mole),  and  the  gas-phase  rate  constant  is  calculated  to 
be  about  ten  times  faster  (1.35  x  10” sec  “*)  using  Rogers'  values*^  (A  * 

3.14  x  10*®  sec-*,  Ea  =  34.1  kcal/mole)  for  the  gas-phase  decomposition. 
However,  this  is  in  agreement  with  Robertson’s  statement’  that  his  manometric 
results  were  Independent  of  sample  weight  and  pressure  of  overlying  gas. 

Thus,  it  appears  that  vapor  reaction  probably  does  not  introduce  as  much  error 
into  non-thermal  kinetic  studies  as  into  thermal  analysis  studies. 


1  1 

R.  N.  Rogers,  private  communication,  Los  Alamos  Scientific  Laboratory,  1979. 

*^R.  N.  Rogers  and  G.  W.  Daub,  "Scanning  Calorimetric  Determination  of  Vapor- 
Phase  Kinetic  Data,"  Anal.  Chem, .  Vol.  45,  p.  596,  1979. 
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Possible  reasons  for  this  discrepancy  include  the  following:  (a)  The  DSC 
experiment  follows  the  decomposition  by  its  heat  release.  If  the  gas-phase 
decomposition  had  a  greater  heat  release  per  mole  than  the  liquid-phase 
decomposition,  the  DSC  apparatus  would  sec  the  gas-phase  reaction  as  being 
proportionately  more  important  than  would  the  conventional  kinetic 
experiments,  thereby  necessitating  a  larger  correction,  (b)  The  m/V  ratio 
for  the  DSC  experiments  (ca  0.5mg/0.022ml*^>*4  23mg/ml  is  generally  larger 
than  those  for  conventional  kinetic  experiments  (see  for  example  Reference  4) 
(ca  40mg/26ml  »  1.5mg/ml);  this  might  promote  any  exothermic  reaction  of 
gaseous  products  with  each  other  or  with  unreacted  RDX,  since  the  higher  ra/v 
ratio  might  tend  to  keep  gaseous  decomposition  products  in  contact  with 
unreacted  liquid  (see  discussion  under  auto-acceleration).  This  effect  might 
conceivably  operate  even  though  gaseous  products  are  generally  allowed  to 
escape  through  a  small  hole  in  the  top  of  the  DSC  pan.  (c)  The  DSC  apparatus 
might  detect  heat  output  from  reaction  of  gaseous  decomposition  products  with 
each  other  rather  than  just  that  from  decomposition  of  RDX  itself. 

The  possibility  should  be  kept  in  mind  that  some  of  the  heat  evolution 

from  the  gas  phase  reaction  may  be  due  to  reaction  among  the  products  rather 

than  just  decomposition  of  RDX;  however,  good  first  order  plots  and  reaction 
orders  of  ^a  0.9  were  obtained  for  the  gas  phase  reaction  so  it  seems  likely 
that  the  observed  rate  constants*4  and  Arrhenius  parameters  (Appendix  A) 

obtained  by  DSC  are  actually  those  for  the  gas  phase  decomposition  of  RDX  and 

HMX,  although  it  is  difficult  to  rule  out  the  possibility  of  pseudo-first- 
order  reaction  of  products.  Expulsion  of  unreacted  HMX  or  RDX  vapor  through 
the  hole  in  the  lid  of  the  pan  might  also  lead  to  artificially  high  rate 
constants. 

The  above  discussion  applies  specifically  to  RDX  but  since  HMX  also 
exhibits  a  vapor  phase  reaction  it  should  apply  to  HMX  also.  However,  the 
situation  may  be  more  complicated  for  HMX;  see  discussion  in  Appendix  B. 

D.  Reactions  on  the  Container  Surface 


In  connection  with  g|s-phase  decomposition,  another  possible  (and  often- 
overlooked)  complication*  is  the  fact  that  reactions  assumed  to  take  place 
entirely  in  the  homogeneous  gas  phase  sometimes  in  fact  take  place  on  or  are 
catalyzed  by  the  wall  of  the  reaction  vessel.  More  work  is  needed,  but  wall 
reactions  are  probably  not  a  very  Important  source  or  error  in  studies  of  gas- 
phase  decomposition  of  H>K  and  RDX  themselves.  This  tentative  conclusion  is 


> 


I  s 

P.  J.  Robinson  and  K.  A.  Holbrook,  "Unimolecular  Reactions,"  Wlley- 
Intersclence.  New  York,  pp.  6-7,  1972. 
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based  on  two  observations:  First,  it  has  been  reported11*  that  the  rate  of 
gas-phase  HMX  pyrolysis  is  independent  of  the  nature  of  the  container  surface 
(note  however,  that  dimethylnitramine  decomposition  did  exhibit  such  a 
dependency^);  second,  consider  the  gas-phase  decomposition  data  for  RDX 
plotted  in  Figure  2.  As  stated  in  Appendix  A  and  in  the  section  on 
"Statistical  Analysis,"  the  data  plotted  includes  DSC  data  of  References  A-l 
and  A-4,  and  manoraetric  data  from  Reference  A-3.  The  DSC  data  were  done  in 
metal-walled  DSC  pans.  The  wall  material  of  the  manometric  (Reference  A-3) 
studies  was  not  stated,  but  was  presumably  some  kind  of  glass  or  quartz.  If 
this  is  the  case,  the  close  correspondence  (Figure  2)  of  data  with  the  two 
types  of  wall  material  would  seem  to  suggest  that  wall  material  doesn't  make 
much  difference,  and  hence  that  wall  reactions  are  not  to  important.  Note 
again,  however,  that  the  SRI  workers**’  found  that  the  rate  of 
dimethylnitramine  decomposition  did  depend  on  wall  material;  thus  wall 
material  is  apparently  important  at  least  for  the  case  of  this  one 
nitramine.  Note  also  the  discrepancy  between  dimethylnitramine  and  HMX  in 
this  respect;  possible  explanations  include  structural  effects  on  absorption 
characteristics,  differences  in  decomposition  mechanism,  and  difference  in 
pressure  at  which  the  experiments  were  carried  out.  See  also  paragraph  6  of 
Note  Added  in  Proof  No.  3. 


This  apparent  importance  of  wall  reactions  for  at  least  some  nitramines 
should  be  kept  in  mind  when  reading  previous  work  on  gas-phase  decomposition 
of  nitramines.  Furthermore,  it  is  recommended  that  future  workers  on 
decomposition  of  nitramines  under  conditions  where  gas-phase  decomposition  may 
be  important,  should  at  least  check  for  wall  reactions  by  examining  the  effect 


of  added  glass  beads,  glass  wool,  etc.,  if  practicable  it  would  probably  also 
be  a  good  idea  to  study  the  effect  of  variations  in  wall  material.1-’  As 


discussed  above,  it  appears  that  HMX  and  RDX  themselves  may  not  be  affected 
too  much  by  wall  reactions,  but  more  information  is  needed  here  also. 


E.  Autoacceleration,  Chain  Reaction  and  Cage  Effects 


In  the  homogenous  liquid  phase,  autoacceleration  of  decomposition  by 
products  might  be  swallowed  up  in  the  overall  first  order  rate  constant  if  the 
autoacceleratory  follow-up  reactions  were  rapid  and  the  chain  lengths 
relatively  short.  Robertson^  reported  some  evidence  of  autocatalysis  of  RDX 
decomposition,  apparently  by  a  stable  product  formed  in  solution. 
Autoacceleratory  behavior  might  in  principle  be  expected  to  either  raise  or 
lower  the  apparent  activation  energy  depending  on  the  relative  temperature 
dependencies  of  the  initial  cleavage  and  autocatalyic  follow-up  steps,  but 
lowering  seems  to  be  the  usual  pattern.  Evidence  for  some  effect  on  liquid 
decomposition  by  gaseous  products  is  provided  by  observations  that  the 
activation  energies  for  RDX  decomposition  were  about  15  kcal/mole  lower  for 


*^D.  F.  McMillen,  J.  R.  Barker,  K.  E.  Lewis,  P.  L.  Travor  and  D.  M.  Golden, 
"Mechanisms  of  Nitramine  Decomposition:  Very-Low  Pressure  Pyrolysis  of  HMX 
and  Dimethylnitramine,"  Final  Report,  SRI  Project  "YU-5787,  18  June  1979 
(SAN  0115/117).  (Supported  by  Department  of  Energy  through  Lawrence 
Livermore  Laboratories). 
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closed  pan  (sealed)  chan  for  open  pan  (l  hole)  DSC  experiments.*^  Note 
however  the  possibility  of  heat  effects  due  to  confinement  of  reactive  gaseous 
products.  In  general,  however,  the  agreement  of  the  most  reliable  rate 
constants  and  activation  energies  with  each  other  and  with  the  vapor-phase- 
reaction-corrected  DSC  results  suggests  that  under  these  conditions  of  less 
confinement  (lower  m/V  ratio),  rate  constants  are  not  being  affected  by 
reaction  of  gaseous  products  with  the  liquid  phase.  It  is  also  possible  that 
case  effects  could  affect  rates  and  activation  parameters.  See  also  Section  V 
A  and  Note  Added  in  Proof  No.  2. 

F.  Isothermal  vs.  Nonisothermal  Thermal  Analysis 

Two  general  approaches  are  generally  used  in  measuring  activation 
energies  and  frequency  factors  by  thermal  analysis  methods.  We  will  refer  to 
these  as  isothermal  and  nonisothermal  approaches.  In  the  first  (isothermal) 
approach  the  decomposing  sample  is  held  at  constant  temperature  and  the  plot 
of  heat  evolved  vs  time  is  used  to  derive  the  rate  constant  at  that 
temperature;  the  rate  constants  derived  in  this  way  at  several  temperatures 
are  then  used  in  a  conventional  Arrhenius  plot  of  log  of  k  vs  l/T  to  derive  Ea 
and  A.  In  the  second  (nonisothermal)  approach  the  activation  energy  is 
derived  by  any  of  a  variety  of  calculation  procedures  from  time-property 
curves  derived  from  programmed  temperature  runs. 

The  isothermal  method  seems  capable  of  giving  good  results  for  the  liquid 
phase  reaction  provided  a  correction  is  made  for  the  gas  phase  reaction  which 
occurs  simultaneously  with  liquid  phase  reaction;  thus  the  activation  energy 
and  frequency  factor  found  for  decomposing  liquid  RDX*^  were  43.1  kcal/mole 
and  2.44  x  10*®  sec"1  respectively  without  corrections  for  vapor  phase 
reaction,  and  47.1  kcal/mole  and  2.02  x  10*®  sec-*  with  correction  for  gas 
phase  reaction,  the  latter  pair  of  values  being  in  reasonable  agreement  with 
Robertson's  values*  of  47.5  kcal/mole  and  3.16  x  10*®  sec-*  (not  2.17  x  10*® 
sec-*  as  stated  in  Reference  12).  As  far  as  this  writer  is  aware,  this  is  the 
only  attempt  at  correcting  isothermal  DSC  data  on  liquid  HMX  and  RDX  for 
concurrent  vapor  phase  reaction;  since  the  correction  in  activation  energy 
seems  to  be  about  47.1  -  43.1  »  4.0  kcal/mole  and  other  isothermal  DSC  values 
(Appendix  C,  References  13,  14,  16,  20,  28-30)  are  generally  in  the  range 
40-44  kcal/mole  it  would  appear  that  application  of  this  correction  would 
raise  these  value 
Robertson's  value 
correction  factor 

whenever  possible  values  of  Ea  and  A  measured  by  thermoanalytical  methods  for 
liquid  phase  HMX  and  RDX  should  be  corrected  for  concurrent  vapor  phase 
decomposition. 

In  the  case  of  the  nonisothermal  approach  it  is  hard  to  know  how  to 
include  a  correction  for  the  gas  phase  reaction.  In  addition  these  methods 
are  generally  subject  to  errors  both  experimentally  and  because  of  assumptions 
Introduced  in  the  course  of  deriving  the  equations  used  in  the  calculation 


to  the  range  44-48  kcal/mole  in  better  agreement  with 
of  47.5  kcal/mole.  A  similar  treatment  suggests  that  the 
in  A  is  about  10*  or  10®;  thus  it  is  recommended  that 


*^K.  Kishore,  "Thermal  Decomposition  Studies  on  Hexahydro-1 ,3 ,5-trinitro-s- 
trlazine  (RDX)  by  Differential  Scanning  Calorimetry,"  Propellants  and 
Explosives.  Vol.  2,  p.  78,  1977. 
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(see  for  example  References  18-21).  Because  of  these  considerations,  and  due 
to  time  limitations  only  a  short  discussion  of  non-isothermal  DSC 
determinations  will  be  given  here.  However,  an  attempt  has  been  made  at 
including  in  the  Appendices  all  available  values  measured  by  these  methods. 

For  example,  it  seems  appropriate  to  point  out  that  it  appears10  that 
one  common  method2 *  for  calculating  activation  energies,  namely  the  method  of 
Kissinger,2*  may  be  subject  to  errors,  when  used  to  study  HMX  and  RDX 
decomposition.  In  particular  Rogers  and  Smith  reported2®  that  the  Kissinger 
method  gave  low  values  for  activation  energy  (28  kcal/mole)  and  frequency 
factors  (log  A(sec~*)  -  10.6)  although  several  o*her  methods  gave  values  more 
in  agreement  with  Robertson.’  However,  Kishore17*22  fou  id  that  the  Kissinger 
method  gave  Ea  ”  42  kcal/mole,  in  approximate  agreement  with  isothermal  and 
other  scanning  methods.  Thus,  it  would  appear  that  until  this  situation  is 
better  understood,  the  Kissinger  method  and  other  nonisothermal  methods  should 
be  used  with  caution;  care  snould  be  taken  that  the  system  under  study 
conforms  to  the  assumptions  made  in  deriving  the  method.  Actually,  the 
physical  coraplexities  alone  may  well  render  HMX  and  RDX  decomposition 
unsuitable  for  study  by  many  nonisothermal  methods. 

It  also  seems  worth  mentioning  that  several  papers17, 20,23  inciU(je 
comparisons  of  Arrhenius  parameters,  obtained  by  various  dynamic  or 
nonisothermal  thermoanalytical  methods,  including  DSC  and  TGA  results.  In 
general,  the  dependence  of  the  results  of  these  studies  on  experimental 
conditions  and  method  of  calculation  appears  to  be  such  that,  as  suggested  by 
Garn2^  (for  systems  other  than  HMX  and  RDX),  it  actually  seems  more 
appropriate  to  regard  them  as  "Temperature  Coefficients  of  Reaction"  rather 
than  “Activation  Energies,"  since  they  almost  certainly  do  not  represent  an 
energy  of  reaction  for  a  single,  homogeneous  chemical  process,  but  rather  the 


1«W.  W.  Vendlandt,  "Thermal  Methods  of  Analysis,"  2nd.  Edition,  John  Wiley  and 
Sons,  New  York,  NY,  p.  187-193,  1974. 

*^J.  H.  Sharp,  "Differential  Thermal  Analysis,"  R.  C.  Mackenzie,  ed.,  Vol.  2, 
Chapter  28,  See  especially  p.  35,  Academic  Press,  New  York,  NY,  1972. 

20R.  N.  Rogers  and  L.  C.  9mith,  "Application  of  Scanning  Calorimetry  to  the 
Study  of  Chemical  Kinetics,”  Thermochlmlca  Acta,  Vol.  I,  p.  1,  1970. 

«H.  E.  Kissinger,  "Reaction  Kinetics  in  Differential  Thermal  Analysis,”  Anal. 
Chem..  Vol.  29,  p.  1702,  1957. 

22K.  Kishore,  "Comparative  Studies  on  the  Decomposition  Behavior  of  Secondary 
Explosives  RDX  and  HMX,"  Defense  Science  Journal,  Vol.  18,  p.  59,  See  Chem. 
Abstr,,  Vol.  90,  206758V,  1978. 

“JB.  D.  Smith,  "Differential  Thermal  Analysis  and  Dynamic  Thermogravimetric 
Analysis  Studies  of  RDX  and  HMX,"  NWL  Technical  Report  TR-2316,  July  1969 
(AD-857  655). 

2^P.  D.  Gam,  "Temperature  Coefficient  of  Reaction,”  Thermochlmlca  Acta, 

Vol.  28,  p.  185,  1977. 


temperature  dependence  of  a  complex  process  Involving  both  vapor  and  condensed 
phase,  and  most  likely  some  degree  of  autocatalysis  as  well.  Thus,  although 
these  methods  may  well  have  usefulness  in  studying  trends  in  data  for  a  series 
of  experiments  under  similar  conditions,  the  absolute  values  probably  have 
meaning  only  under  the  exact  conditions  for  which  they  were  determined. 

III.  STATISTICAL  ANALYSIS 

The  statistical  analysis  was  carried  out  by  two  methods.  The  first 
method  was  the  method  of  Cvetanovic  and  Singleton^ »26  which  involves  fitting 
the  logarithms  of  the  data  using  a  transformation  of  the  weights.  The  result 
•  a  least  squares  fit  of  the  original  data.  The  program  used  had  been 
written  for  another  project  by  Dr.  T.  P.  Coffee,  and  a  copy  of  it  is  included 
as  Appendix  E.  The  second  program  fits  data  to  the  logarithmic  form  (log  kc  * 
log  A  -  Ea/RT)  of  the  Arrhenius  Equation.  The  result  is  a  least  squares  fit 
of  the  logarithms  of  the  data.  This  program  was  written  by  D.R.  Crosley  and 
G.E.  Keller  for  other  work.  A  copy  of  this  program  is  included  as  Appendix  F 
of  the  present  report.  The  author  thanks  Drs.  Coffee  and  Keller  for  carrying 
out  the  analyses  of  the  data  using  their  respective  programs. 

Data  on  gas-phase  decomposition  of  HMX  and  RDX  was  taken  from  Tables  A-2 
and  A-3  of  Appendix  A.  These  Tables  include  all  data  known  to  the  present 
writer  on  gas-phase  decomposition  of  HMX  and  RDX.  Data  analyzed  for 
decomposition  of  neat  liquid  HMX  were  taken  from  Table  &-3  of  Appendix  B; 
these  are  believed  by  the  present  writer  to  be  the  most  reliable  data  on 
decomposition  of  neat  liquid  HMX.  Although  all  three  sets  of  data  were 
analyzed  and  plotted  together,  it  is  suggested  that  until  the  question  of  the 
gas-phase  correction  for  HMX  decomposition  is  better  understood,  the  best 
value  for  liquid  HMX  would  probably  be  those  derived  from  Robertson's  work 
(Reference  9).  The  data  analyzed  for  decomposition  of  neat  liquid  RDX  were 
taken  from  Table  C-3;  for  reasons  discussed  in  Appendix  C,  these  are  believed 
by  the  present  writer  to  be  the  best  data  for  decomposition  of  neat  liquid 
RDX.  Since  several  authors  reported  data  (Table  C-4)  on  decomposition  of  RDX 
in  TNT  solution,  these  data  were  analyzed  and  plotted  together  also.  The  data 
for  neat,  liquid  RDX,  gaseous  RDX,  neat  liquid  HMX,  gaseous  HMX,  and  RDX  in 
TNT  solution  are  plotted  in  Figures  1-5  respectively,  and  the  Arrhenius 
parameters  obtained  from  the  two  programs  are  summarized  in  the  first  four 
columns  of  data  in  the  Table.  The  last  two  columns  of  data  in  the  Table  list 
the  values  suggested  by  the  present  writer  for  decomposition  of  the  indicated 
compounds  in  the  indicated  state  of  aggregation.  In  general,  these  are  the 
values  from  the  Cvetanovic-Singleton  Program,  as  it  seems  that  the  direct  fit 
to  the  logarithmic  form  should  be  more  accurate. 25.26 


7  r 

R.  J.  Cvetanovic  and  D.  L.  Singleton,  "Comment  on  the  Evaluation  of  the 
Arrhenius  Parameters  by  the  Least  Squares  Method,"  Int.  J.  Chem.  Kinet., 
Vol.  9,  p.  481,  1977. 

26 

R.  J.  Cvetanovic  and  D.  L.  Singleton,  "Comment  on  the  Evaluation  of  the 
Arrhenius  Parameters  by  the  Least  Squares  Method,"  Int.  J.  Chem.  Kinet.. 
Vol.  9,  p.  1007,  1977. 
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TABLE  1.  SUGGESTED  LOGlnA(sec-1)  AND  E 
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Plot  of  Log  k  vs  1000/T  or  Liquid  RDX  Decomposition. 
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Figure  2.  Plot  of  Log  k  vs  1000/T  for  Decomposition  of  Gaseous  RDX. 
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Figure  4.  Plot  of  Log  k  vs  1000/T  for  Gas-Phase  HMX  Decomposition 


Figure  5.  Plot  of  Log  k  vs  1000/T  for  Decomposition  of  RDX 

in  Solution  in  TNT. 


The  values  of  logjo  A  listed  in  the  Table  are  given  to  three  significant 
figures  past  the  decimal  point,  although  it  is  highly  questionable  whether  the 
accuracy  of  the  intercepts  themselves  justifies  more  than  one  such  figure. 

The  other  figures  are  given  because  it  is  felt  that  modelers  may  find  more 
precision  useful  in  interpolating  within  the  temperature  range  covered  by  the 
actual  measurements,  or  for  extrapolating  to  temperatures  not  too  far  outside 
these  ranges.  Note  that  the  original  rate  constants  and  temperatures  on  which 
the  values  given  arc  based,  are  for  the  most  part  given  to  three  or  four 
significant  figures  by  the  original  authors. 


IV.  GENERAL  DISCUSSION  AND  SUGGESTED  VALUES  FOR  KINETIC  PARAMETERS 
A.  Gas-Phase 

There  is  a  need  for  more  work  on  the  gas-phase  kinetic  parameters  for  HMX 
and  RDX  decompositions,  as  mentioned  in  the  discussion  of  data  in  Appendix  A. 
The  best  values  for  use  in  modeling  studies  involving  gas-phase  decomposition 
at  relatively  low  temperatures  v.ca  200-300°C)  would  probably  be  log  A  (sec-*) 

=  11.991  and  Ea  =  30.5  kcal/mole  for  RDX,  based  on  the  plot  and  statistical 
analysis  of  all  data  (Table  and  Figure  2).  For  HMX,  the  data  is  too  scattered 
to  allow  choice  of  a  single  recommended  value  (Table,  Figure  A)  but  the  values 
are  most  likely  in  the  range  30-35  kcal/mole  and  log  A  (Sec-*)-  «  12.8 
respectively.  Until  more  data  is  available,  Ea  *  32.5  kcal/mole  would 
probably  be  as  good  as  anything. 

Furthermore,  there  is  reason  to  suspect  that  as  temperature  rises  the 
Arrhenius  parameters  may  change  appreciably,  due  to  a  change  in  decomposition 
n.jchanism.  This  follows  from  the  report  by  McMillen  et  al.,*^  who  carried  out 
thermochemical  estimations,  very  low  pressure  pyrolysis  and  mas:*  spe-cral 
studies  on  gas  phase  HMX  and  dimethyi-N-nitroamine  decomposition;  their 
results  were  in  agreement  with  the  idea  that  at  lower  temperatures  the 
principal  reaction  in  the  vapor  phase  was  HONO  elimination,  and  at  higher 
temperatures  N-N02  cleavage,  with  the  cross-over  point  occurring  at  about 
500°K  or  300°C.  If  this  is  the  case,  it  might  be  a  good  Idea  if  possible  to 
represent  the  low  temperature  reaction  by  the  values  given  above.  The  high 
temperature  N-N02  cleavage  reaction  could  be  represented  by  rhe  parameters 
estimated  by  Shaw  and  Walker^'  for  N-N02  cleavage  of  HMX:  log  A  (sec-*)  * 

16.4  and  Ea  »  46.2  kcal/mole.  For  RDX  Ea  should  be  about  the  same  but  log  A 
might  be  expected  to  be  lowered  from  16.40  to  16.28  since  RDX  has  only  three 
equivalent  N-N02  groupings.  If  possible  it  would  probably  be  best  for 
modelers  to  incorporate  this  possible  mechanism  shift  by  representing  low 
temperature  reactions  as  proceeding  on  the  basis  of  whatever  temperature ~is 
predicted  by  the  model  under  consideration.  Note  also  that  the  exact  location 
of  the  crossover  point  is  uncertain;  thus  Shaw  and  Walker,  '  using  a  somewhat 
different  set  of  parameters,  got  a  crossover  point  of.  ca  60°C.  In  closing,  we 


2^R.  Shaw  and  F.  E.  Walker,  "Estimated  Kinetic  and  Thermochemistry  of  Some 
Initial  Unimolecular  Reactions  in  the  Thermal  Decomposition  of  1,3,5,7- 
Tetranitro-1 ,3,5,7-tetraazacyclooctane  in  the  Gas  Phase,"  J.  Phys.  Chem. 
Vol.  81,  p.  2572,  1977. 
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point  out  again  that  there  is  a  need  for  more  data  on  the  gas-phase 
decomposition  of  both  HMX  and  RDX,  especially  in  the  high-temperature  range, 
where  N-NO2  cleavage  may  become  important. 

Available  data  on  gas-phase  decomposition  of  nitramines  other  than  HMX 
and  RDX  has  been  summarized  by  the  SRI  group, ^  who  also  carried  out 
additional  measurements  themselves,  on  dimethylnitramine  and  HMX.  The 
literature^--^  activation  energies  and  frequency  factors  varied  widely  (Ea 
values  reported  for  dimethylnitramine  varied  from  37  to  53  kcal/mole,  for 
example),  and  were  generally  interpreted  as  being  consistent  with  log  A(sec“*) 
■  12.4  and  Ea  ■  35-39  kcal/mole,  or  with  log  A(sec~*)  »  16.5  and  Ea  *  43-48 
kcal/mole.  These  values  were  of  interest  because  they  are  characteristic  of 
HONO  elimination,  and  of  N-NO2  cleavage  respectively.  The  SRI  workers 
interpreted  their  results  in  terms  of  temperature-variation  of  the  relative 
importance  of  these  two  mechanisms,  see  the  preceding  paragraph.  However,  tne 
role  of  wall  reactions  seems  uncertain  at  this  point.  As  pointed  out  by  the 
SRI  workers, ^  the  apparent  lack  of  surface  effects  on  HMX  decomposition  seems 
puzzling,  since  dimethylnitramine  decomposed  faster  with  a  normally  aged  or 
oxygen  cleaned  surface  than  with  a  carbon-coated  surface.  On  the  other  hand, 
the  decomposition  rate  of  N-methyl-N-choloromethylnitramine  was  unaffected  by 
increasing  the  specific  surface  area  of  the  container  by  over  an  order  of 
magnitude, and  Figure  1  of  Reference  30  suggests  that  decomposition  kinetics 
of  dlethylnltramine  were  not  greatly  affectea  by  filling  the  reaction  vessel 
w..n  glass  packing.  The  explanation  for  these  results  seems  uncertain; 
differences  in  pressure  between  experiments  may  make  a  difference,  or 
pos8iblly  the  surface  reaction  has  a  sterlc  requirement  that  is  satisfied  only 
by  the  relatively  unhindered  dimethylnitramine.  In  any  case,  there  Is  clearly 
a  need  for  further  work  on  the  nature  and  occurrence  of  surface  reactions  in 
the  gas  phase  decomposition  of  gaseous  nitramines,  if  the  gas-phase 
decomposition  Itself  is  to  be  understood. 

B.  Liquid  Phase 


The  activation  energy  and  frequency  factors  from  the  four  apparently  most 
reliable  data  sources  on  decomposition  of  liquid  RDX  have  been  described  in 
Appendix  C  and  accompanying  discussion.  A  plot  of  the  data  for  all  four 
sources  on  one  axis  is  given  in  Figure  1,  and  a  statistical  analysis  of  the 


M.  Fluornoy,  "Thermal  Decomposition  of  Caseous  Dimethylnitramine,"  J. 
Chem.  Phys..  Vol.  36,  p.  1106,  1962. 
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B.  L.  Korsunskii,  and  F.  I.  Dubovitskii,  "Thermal  Decomposition  Kinetics  of 
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J  B.  L.  Korsunskii,  F.  I.  IXibovitskii ,  and  E.  A.  Shurygin,  "Kinetics  of  the 
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points  for  all  four  sotirces  by  the  procedure  of  Cvetanovlc  and  Sineleton 
(pp.  19-26)  gave  A  =  2.46  x  10^  (log  A  =  18.668)  and  =  47.8  kcal/mole: 
these  are  our  suggested  values  for  decomposition  of  neat  liquid  RPX. 


The  values  for  neat  liquid  HMX  were  similarily  discussed  in  the 
discussion  accompanying  Appendix  B.  The  most  reliable  values are  those 
chat  were  plotted  (Figure  3),  and  analyzed  similarily  to  the  FPX  values  giving 
values  for  log  A  and  F,^  which  could  possibly  be  presented  as  recommended 
values  for  modeling  studies  involving  decomposition  of  pure  liquid  HMX.  Mote, 


to  appear 
this 

situation  is  understood  better  Jjhan  at  present,  it  may  be  a  good  idea  to  just 
use  Robertson's  manometric  data5  as  the  best  data  for  decomposition  of  neat 
liquid  HMX:  as  shown  in  Appendix  B.  this  gives  E  =  53.1  kcal/mole  and  log  A 
(sec  =  19.950.  Note  also  that  the  values  of  log  A  (sec  )  (=  21-21.5) 
obtained  from  all  three  sets  of  data  together  seem  unreasonably  high. 
Alternatively,  Robertson's  values  of  =  52.7  kcal/mole  and  log  A  (sec  )  = 
19.7  could  be  used,  especially  since  some  of  the  discrepancy  between  the  two 
sets  was  presumably  introduced  by  the  present  writer  in  the  course  of  reading 
numbers  off  Robertson's'  published  graph. 


32r.  J.  Powers,  AFATL,  Eglin  AFB ,  FL,  Private  Communication,  1980. 
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Activation  energies  for  thermal  decomposition  of  nitramines  other  than 
HMX  and  RDX^»34-39  t^e  neat  state  are  generally  in  the  range  40±5 

kcal/mole,  with  frequency  factors  correspondingly  decreased.  The  values  have 
been  observed  to  decrease  with  increasing  confinement  in  at  least  one  case,  5 
namely  that  of  di-( 2-nitroxyethyl)-N-nitroamine;  this  was  attributed^  to 
autoacceleration  due  to  a  combination  of  hydrolysis  ard  oxidation  involving 
water  and  NO2  formed  in  the  decomposition. 

C.  Solution-Phase  Decomposition  of  RDX  and  HMX 

Activation  energies  and  log  A  for  solution  phase  decomposition  of  HMX  and 
RDX  are  summarized  in  Appendix  B  (HMX)  and  Appendix  C  (RDX).  Like  the  gas- 
phase  values  the  values  for  frequency  factors  and  activation  energy  are 
generally  lower  for  decomposition  in  solution  than  in  pure  liquid.  There  was 
enough  data  on  the  decomposition  of  RDX  in  solution  in  TNT  to  allow 
statistical  analysis  by  the  procedure  of  Cvetanovic  and  Singleton;  a  plot  is 
given  in  Figure  5.  The  line  for  the  three  studies  at  concentration  of  20-60% 
falls  above  the  line  defined  by  the  studies  of  Robertson  (carried  out  at 
concentrations  of  1-5%  RDX)  within  the  region  studied  experimentally,  although 
the  frequency  factor  for  the  20-60%  points  is  lower  than  that  for  the  1-5% 
points;  apparently  there  is  a  line  crossing  resulting  from  a  slightly  higher 
activation  energy  for  the  1-5%  studies.  It  is  not  certain  whether  tnis  line 
crossing  is  real  or  is  an  artifact  resulting  from  scatter  in  the  experimental 
data.  If,  as  it  does  not  seem  altogether  unreasonable,  the  line-crossing  is  a 
artifact,  then  it  would  appear  that  the  studies  at  20-60%  concentration  have  a 
frequency  factor  higher  by  a  factor  of  ca  2-3  than  the  studies  at  1-5%  RDX. 


-*4F.  I.  Dubovitskii,  0.  B.  Manelis  and  L.  P.  Smirnov,  "Kinetics  of  the  Thermal 
Decomposition  of  N-Methyl-N,  2,  4,  6-Tetranitroaniline  (Tetryl),"  Russian 
Journal  of  Physical  Chemistry.  Vol.  35,  p.  255,  1961. 

„  ^5p.  I.  Dubovitskii,  Yu.  I.  Rubtsov,  V.  V.  Barykin  and  G.  B.  Manelis, 

|  “Kinetics  of  Thermal  Decomposition  of  Diethylnitramine  Dinitrate,"  Bull. 

Acad.  Sci.  USSR,  Div.  Chem.  Sci.,  p.  1126,  1960. 

^B.  S.  Svetlov  and  B.  A.  Lur'e,  "Thermal  Decomposition  of  Di(nitroxyethyl) 
nitramine,"  Russian  Journal  of  Physical  Chemistry,  Vol.  37,  p.  1073,  1963. 
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However,  the  possibility  of  a  systematic  error  in  the  1-5%  studies  should  be 
kept  in  mind,  since  these  studies  are  all  taken  from  a  single  publication. 
Continuing  this  possible  trend,  the  frequency  factor  for  decomposition  of  neat 
liquid  RDX  is  considerably  higher  than  for  the  reactions  in  solution  (2.46  x 
10f8  sec-*,  Figure  1).  The  values  for  gas-phase  decomposition  (Table,  p.  26- 
27)  are  also  lower  than  tor  the  neat  liquid;  possible  reasons  for  this 
increase  in  frequency  factor  with  concentration  include  (a)  a  difference  in 
mechanism  of  energy  transfer  leading  to  an  increase  with  RDX  concentration  in 
the  efficiency  of  energy  transfer  into  RDX  molecules;  (b)  an  increase  with  RDX 
concentration  in  the  importance  of  biraolecular  follow-up  steps  much  as  those 
discussed  in  our  discussions23  of  possible  mechanisms  involved  in  HMX  and  RDX 
decomposition;  and  (c)  catalysis  of  TNT  decomposition  by  solute  RDX.  0 

For  both  HMX  and  RDX,  the  activation  energies  for  decomposition  in 
solution  (Appendices  B,  C  and  Figure  5)  are  generally  lower  (ca  40-45 
kcal/mole)  than  those  (ca  50  kcal/mole)  for  decomposition  in  the  neat  liquid 
phase  (Appendices  B,  C  and  Figures  1  and  3).  The  best  values  for  the  vapor- 
phase  decomposition  of  both  compounds  are  also  lower  than  for  the  neat 
liquid.  Possible  reasons  for  this  include  the  following:  (a)  the  mechanism 
changes  from  one  unimolecular  mechanism  to  another  on  going  from  gas  phase  or 
solution  to  neat  liquid  HMX;  (b)  the  mechanism  stays  the  same,  but  neat  liquid 
HMX  and  RDX  exerts  a  different  solvent  effect  than  the  inert  solvents  studied; 
(c)  the  mechanism  shifts  from  unimolecular  in  the  gas-phase  or  in  dilute 
solutions,  to  bimolecular,  possibly  chain  or  autocatalytic,  in  neat  liquid  HMX 
or  RDX.  Possible  bimolecular  steps  might  include  reaction  of  a  nitro  group  on 
one  molecule  of  HMX  or  RDX  with  a  CH2  grouping  of  another  HMX  or  RDX  molecule, 
as  illustrated  in  Scheme  VI  of  Reference  2-A,  or  electron  transfer  between  two 
molecules  of  HMX  or  RDX.  A  bimolecular  first  step  for  the  decomposition  in 
the  neat  liquid  is  not  Inconsistent  with  the  observed  first-order  kinetics, 
since  the  concentration  of  HMX/RDX  seen  by  any  given  molecule  of  HMX/RDX  would 
remain  constant  throughout  the  reaction;  hence  the  reaction  would  be  zero- 
order  in  the  second  molecule  of  HMX/RDX  and  first-order  overall.  However,  (c) 
seems  unlikely  in  view  of  the  high  (log  A  -  1018  -  1020)  frequency  factors  for 
the  liquid-phase  reaction,  and  of  the  apparent  first-order  kinetics  in 
solution.  See  also  the  discussion  of  cage  effects  in  Note  Added  in  Proof  No. 
2. 

In  connection  with  the  discussion  of  the  decomposition  of  RDX  in  TNT 
solution,  it  seems  worth  mentioning  that  while  RDX  catalyzed  the  thermal 
decomposition  °f0TNT,^  the  TNT  acts  as  an  inert  solvent  with  regard  to  the 
decomposing  RDX. 


P.  Colraan  and  F.  C.  Rauch,  "Studies  on  Composition  B,"  Final  Report, 
Contract  No.  DAAA21-70-0531 ,  American  Cyanamid  Company,  February  1971  (AD- 
881  190). 
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The  above  values  for  log  A  and  Efl  are  approximately  in  agreement  with 
those  measured41  for  a  series  of  substituted  dialkylnitroamines  in 
dibutylphthalate  solutions;  these  measurements  gave  log  A(sec-1)  *  13.7-14.4 
and  Ea  "  40-42  kcal/mole.  However  studies  on  a  series  of  4-substituted-? ,6- 
dinitrophenylmethyl-N-nitroamines  in  dinitrobenzene  solution  gave  somewhat 
lower  values:  log  A(sec_1)  *  12.4-13.5  and  Efl  -  32-36  kcal/mole. 

D.  Decomposition  of  ,,Solld“  HMX  and  RDX 

The  reason  for  putting  quotation  marks  around  the  work  "solid"  is  that  in 
most  low  temperature  (100-300®C)  work  the  decomposition  of  FDX  and  (possibly 
with  somewhat  less  certainty)  HMX  below  their  melting  point  may  well  involve 
the  liquid  and  vapor  phases  in  addition  to,  if  not  instead  of,  true  solid- 
state  decomposition.  For  example,  Cosgrove  and  Owen4**-  found  that  the  rate 
of  decomposition  of  RDX  at  195°C  (ca  5-10°  below  its  melting  point)  was  (a) 
directly  proportional  to  the  volume  of  the  reaction  vessel;  (b)  for  a  constant 
volume  independent  of  the  amount  of  RDX  present:  and  (c)  retarded  by  the 
pressure  of  inert  gases,  for  example,  nitrogen.  They  concluded  that  RDX  does 
not  decompose  in  the  solid  state  to  any  significant  degree  at  that 
temperature.  Rather,  in  the  early  stages  of  the  decomposition,  the  FDX  was 
considered  to  vaporize  and  the  vapor  to  decompose  to  liquid  products  which 
dissolve  the  solid  and  thus  accelerate  the  decomposition. 

In  agreement  with  this,  studies^  on  decomposition  of  FDX  at  high  (up  to 
5  kbar)  pressures  indicate  that  decomposition  of  RDX  is  much  slower  in  the 
solid  state  just  below  its  melting  point,  than  in  the  liquid  state  at  its 
melting  point. 

Since  HMX  and  RDX  are  very  similar  chemically,  these  results  raise  strong 
suspicions  that  HMX  may  behave  in  a  similar  manner  when  decomposed  below  its 
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melting  point.  In  fact,  Maksimov2***  found  that  HMX  did  behave  similarly  to  RDX 
when  it  was  .'ecomposed  below  its  melting  point,  although  in  the  case  of  HMX, 
the  volume  dependence  was  not  as  great  as  for  RDX.  Thus  the  true  solid-state 
decomposition  may  well  be  more  important  for  HMX  than  for  RDX,  possibly 
because  HMX  is  solid  at  temperatures  higher  than  the  melting  point  of  RDX,  or 
because  of  lower  vapor  pressure  for  HMX  than  for  RDX  at  the  temperatures  of 
the  respective  experiments. 

However,  there  is  another  point  of  view  in  the  question  of  vapor-phase 
involvement  in  "solid”  HMX/ RDX  decomposition:  Batten^^-^®  has  also  studied 
the  decomposition  of  RDX  below  the  melting  point  and  has  interpreted  his 
results  in  terms  of  positive  or  negative  catalysis  of  the  condensed  phase 
decomposition  by  gaseous  decomposition  products. 

After  the  publication  of  Batten's  work,  Cosgrove  and  Owen*^  argued  again 
in  favor  of  the  gas-phase  decomposition,  saying  that  the  results  quoted  by 
Batten  are  not  at  variance  with  the  suggestion2*^ *43  that  the  gas-phase 
decomposition  is  Important  in  the  early  stages  of  the  reaction.  In  addition 
to  the  arguments  of  Cosgrove  and  Owen,  this  writer  feels  that  It  is  worth 
noting  that  Batten2*'  reported  that  with  0.4  g  of  RDX  in  a  standard  sample 
tube,  the  induction  rate  (in  percent  per  minute)  was  one-half  that  with  only 
0.2  g  RDX  in  the  sample  tube.  Since  the  sample  size  was  twice  as  great,  this 
corresponds  to  a  constant  rate  of  decomposition  in  grams/minute.  Since  the 
same  constant-volume  reactor  was  presumably  used  for  both  experiments,  this 
result  seems  more  consistent  with  vaporization  and  gasj-phase  decomposition  if 
it  is  considered  that  the  second  0.2  g  of  RDX  must  He  underneath  the  first 
0.2  g,  and  hence  its  vaporization  is  suppressed.  In  general,  it  seems  to  this 
writer  that  the  increase  in  rate  with  increasing  reactor  volume,  or  amount  of 
free  space,  spreading,  etc.,  and  its  independence  on  sample  weight  at  constant 
volume,  suggest  that  the  gas  phase  decomposition  is  Important  in  the  early 
stages  of  the  decomposition  below  the  melting  point.  Note  however  that  the 
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decomposition  residue  is  apparently  a  positive  catalyst^ »-* ^  and  that  such 
apparent  products  as  formaldehyde,  NO2  and  hydroxymethylformamide  can  catalyze 
the  reaction  either  positively  or  negatively  depending  on  the  exact  reaction 
conditions  (see  for  example  References  44  and  49).  Thus,  it  may  well  be 
possible  to  explain  just  about  any  set  of  results  by  some  combination  of 
inhibition  and  catalysis  of  the  condensed  phase  reaction  by  gaseous  products 
such  as  formaldehyde  and  NO2 ,  expecially  when  the  possibility  of  reaction 
between  these  products  is  considered. 

Clearly,  there  is  a  need  for  more  work  on  the  relative  roles  of  vapor- 
phase,  liquid-phase,  and  true  solid-state  reaction  in  the  decomposition  of  HMX 
and  RDX  below  their  melting  points,  in  order  that  the  contribution  due  to  true 
solid-state  decomposition  can  be  extracted  and,  extrapolated  Intelligently  to 
combustion  and  explosion  conditions.  One  possible  type  of  experiment  that 
might  help  here  would  be  to  use  various  amounts  of  sample  in  varying  amounts 
of  sample  tubes  of  similiar  construction,  in  reactors  of  differing  volumes  but 
otherwise  similar  construction.  This  would  make  it  possible  to  get  a  better 
idea  of  the  effect  of  sample  size,  configuration,  and  reactor  volume,  each 
with  the  others  held  constant.  This  should  be  done  for  HMX  as  well  as  for 
RDX. 


Thus,  these  below-melting  point  activation  parameters  are  actually 
overall  numbers  for  complex  processes,  and  it  is  difficult  to  relate  them  to 
any  single,  homogeneous  chemical  reaction.  Garn^  has  suggested  the  use  of 
the  term  "Temperature  Coefficient  of  Reaction"  rather  than  "Activation  Energy" 
in  situations  of  this  type;  the  present  writer  agrees  with  this  suggestion. 

Furthermore,  Arrhenius  parameters  for  decomposition  below  the  melting 
point  seem  to  be  dependent  on  such  factors  as  degree  of  spreading  of  the 
sample,  which  are  difficult  to  evaluate  in  view  of  the  limited  infomation 
given  in  many  of  the  reports.  For  example,  Batten  and  Murdie^®  found  that 
activation  energies  for  decomposition  of  RDX  below  its  melting  point  were 
dependent  on  the  degree  of  spreading  of  the  sample. 

For  these  reasons  and  because  of  time  limitations,  a  complete  discussion 
and  evaluation  of  Arrhenius  parameters  for  decomposition  of  nominally  "solid” 
HMX  and  RDX  below  the  melting  point  does  not,  at  this  time,  seem  justified  in 
terras  of  its  relevance  to  any  true  solid-state  decomposition  of  propellants 
that  may  occur  under  combustion  conditions.  Consequently,  although  for  future 
reference  an  attempt  has  been  made  at  including  all  available  activation 
energies  for  decomposition  of  HMX  and  RDX  below  their  melting  points  in 
Appendix  D,  they  will  not  be  discussed  in  detail;  only  a  few  trends  and 
particularly  relevant  points  will  be  considered. 

It  seems  worthwhile  to  begin  the  discussion  of  these  points  by  saying 
something  about  the  relative  rates  of  the  gas,  liquid  and  solid  decomposition 
and  their  relations  to  combustion  modeling  efforts.  The  above-mentioned  work 
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of  Cosgrove  and  Owen^“44  indicates  that  the  gas-phase  reaction  is  much  faster 
than  the  solid-state  decomposition.  Furthermore,  the  liquid-phase 
decomposition  is  faster  than  the  solid-phase  reaction, ^ ^ even  up  to 
pressures  of  5  kbar.^  ^  liquid/solid  rate  ratio  of  10/1  or  greater  seems 
reasonable  for  RDX  at  its  melting  point,  based  on  examination  of  the  pressure¬ 
time  data  in  Figure  8  of  Reference  3.  Also,  decomposition  of  HMX  accelerates 
sharply  when  it  is  heated  past  its  melting  point.  In  view  of  these  results, 
it  seems  quite  reasonable  for  combustion  modelers  to  assume  that  the  solid- 
phase  decomposes  relatively  slowly,  and  that  primary  decomposition  takes 
place,  almost  entirely  in  the  liquid  layer  or  in  the  vapor  phase.  However 
some,  albeit  slow,  solid  state  reaction  apparently  does  take  place, ^  even  for 
RDX  and  pressures  and  temperatures  of  combustion  are  much  higher  than  those  at 
which  the  decomposition  studies  are  carried  out;  it  is  quite  possible  that 
these  higher  temperatures  and  pressures  may  lead  to  accelerated  solid  state 
decomposition.  This  point  seems  especially  relevant  to  conditions  where  there 
is  no  liquid  layer.  It  would  also  help  to  know  whether  at  the  heating  rates 
involved  in  combustion,  HMX  has  time  to  transform  from  the  6-  to  the 
6-polymorph. 


Another  interesting  trend  is  that,  in  many  cases  (Appendix  D)  HMX 
exhibits  a  change  in  activation  energy  in  the  region  240-260°C.  Sometimes 
there  are  two  changes  separated  by  ca  10-20°.  This  break  appears  in  the 
results  of  studies  by  thermal  analysis  and  or  manoraetric  studies,  so  it  is  not 
a  peculiarity  of  any  one  experimental  method.  One  possible  explanation  for 
this  might  be  as  follows:  Goshgarian  (p.  13  of  Reference  11)  has  observed  an 
endotherm  in  DSC  curves  of  HMX  samples,  as  well  as  increased  gas  evolution,  at 
ca  250°C.  He  tentatively  attributes  this  to  melting  of  3-HMX  left  over  from 
the  8-6  phase  transition  at  ^a_  190°C,  apparently  following  an  earlier  report 
by  Tettsov  and  McCrone-^  that  some  8-HMX  survives  the  8-6  phase  transition  anci 
melts  at  240°C;  these  observations  are  described  more  fully  in  Reference  54  of 
the  present  report  than  in  Reference  6  of  Reference  11.  Photographs  of  this 
effect  are  shown  on  pages  19-21  of  Reference  1.  If  the  8-HMX  melts,  it 
presumably  exhibits  a  tendency  to  solidify  into  6-HMX,  seeded  by  the  5-HMX 
already  present.  While  liquid  HMX  was  present  it  might  decompose  at  a  faster 
(preceding  paragraph)  rate  than  the  solid  HMX  and  changes  with  temperature  in 
the  amount  of  liquid  HMX  present  might  cause  a  shift  in  temperature-dependence 
of  the  decomposition  rate,  as  observed.  Even  if  the  8-HMX  melt  hypothesis 
proves  incorrect,  the  activation  energy  shift  might  still  be  related  to  the 
endotherm  and  increase  in  gas  evolution  at  250°C. 


c  2 

T.  B.  Joyner,  "Thermal  Decomposition  of  Explosives.  Part  I.  Effect  of 
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p.  1789,  1977. 

^A.  S.  Teetsov  and  W.  C.  McCrone,  Jr.,  "The  Microscopical  Study  of  Polymorph 
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Another  Interesting  aspect  is  that  BSC  and  Flow  Reactor  Mass 
Spectrometric  (FRMS)  studies  on  pure  HMX  and  RDX  decomposing  at  or  very  close 
to  their  melting  points  (Appendix  D)  often  show  very  high  activation  energies 
of  100-200  kcal/mQle  or  more.  These  have  been  attributed55 ■ 5^  to  thermal 
chain  reactions,  '  or  to  the  complexity  of  the  decomposition  reaction5  or  to 
the  ocurrence  of  processes  involving  the  simultaneous  breaking  of  more  than 
one  chemical  bond,  although  this  last  seems  more  likely  to  lower  than  to 
raise  the  activation  energy.  In  view  of  the  faster  decomposition  rate  for  the 
liquid  than  for  the  solid  state,  the  following  explanation  might  also  be  worth 
considering:  over  the  melting  range  of  HMX  (or  RDX)  a  small  increase  in 
temperature  will  cause  some  of  the  solid  to  melt.  Because  of  the  greater 
decomposition  rate  of  the  liquid,  this  causes  the  rate  to  increase,  resulting 
in  an  artificially  large  temperature-dependence  of  the  decomposition  rate, 
until  the  solid  is  entirely  melted  and  the  temperature  dependence  of  the 
decomposition  becomes  that  of  the  neat  liquid. 


V.  CHEMICAL  MECHANISMS 


Autocatalysis 


When  decomposition  takes  pi ?ce  under  very  confined  conditions,  activation 
energies  for  liquid  decomposition  of  HMX  and  RBX  can  apparently  be  lowered  due 
to  autocatalysis  by  gas-phase  products.  For  example,  when  RDX  was 
decomposed* '  in  an  open  (one  hole  in  lid)  DSC  pan,  F.  was  ^a^  41  kcal/mole  but 
when  the  decomposition  took  place  in  a  closed  pan,  the  rate  was  faster  and 
activation  energy  was  _C£  10-15  kcal/mole  lower.  A  similar  pattern  has  also 
been  observed  for  dinitroxydiethylnitramine;  it  was  attributed  by  the  authors 
to  a  combination  of  hydrolysis  and  oxidation  by  NO2  formed  in  the 
decomposition.  A  similar  effect  may  hold  a  high  pressure,  where  products 
would  have  a  more  difficult  time  escaping  from  the  surface.  Thus,  serious 
consideration  should  be  given  to  the  possibility  that  apparent  activation 
energy  (or  temperature  coefficient  of  reaction)  might  be  lower  under  the 
higher  pressure  characteristic  of  combustion  conditions  than  the  values  given 
above.  However,  it  is  difficult  to  rule  out  the  possibility  that  the  above 
effects  might  be  heat  effects  due  to  confinement  of  reactive  products  near  the 
starting  material.  Also,  under  the  higher  temperatures  of  combustion 
conditions,  higher  activation-energy  unimolecular  reactions  might  predominate. 


5R.  N.  Rogers  and  E.  D.  Morris,  Jr.,  "On  Estimating  Activation  Energies  with 
a  Differential  Scanning  Calorimeter,"  Anal.  Chem. .  Vol.  38,  p.  412,  1966. 

^P.  G.  Hall,  "Thermal  Decomposition  and  Phase  Transitions  in  Solid 
Nitramines,"  J.  Chem.  Soc.,  Faraday  Trans.  Pt.  2,  Vol.  67,  p.  556,  1971. 

^B.  V.  Novazhilov,  "The  Temperature  Dependence  of  the  Kinetic  Characteristics 
of  Exothermic  Reactions  in  the  Condensed  Phase,”  Dokl.  Akad.  Nauk.  SSSR, 
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Dimethylnit ranine  decomposition  also  appeared  subject  to  autocatalysis  by 
Thus,  other  nitramines  also  seem  to  show  autocatalysis. 

Furthermore,  HMX  and  RDX  decompositions  have  been  reported  to  be  accelerated 
or  retarded  by  addition  of  known  decomposition  products  such  as  formaldehyde 
and  oxides  of  nitrogen. ^ 

There  are  also  a  number  of  reports  (see  for  example  References  9,  11,  49, 
51)  of  residues  or  stable  products  that  are  capable  of  accelerating  the 
decomposition  or  of  lowering  the  decomposition  temperature. 

More  work  is  needed  in  the  area  of  autocatalysis  of  nitramine,  HMX  and 
RDX  decomposition.  Since  may  of  the  references  cited  in  the  proceeding 
paragraph  do  not  specifically  allude  to  the  possibility  of  surface  catalysis 
of  gas-phase  reactions,  and  since  the  SRI  workers16  found  that  at  least  some 
of  the  rates  were  affected  by  the  nature  of  the  container  surface,  it  would  be 
helpful  to  have  a  better  understanding  of  the  role  of  surface  reactions  here. 


B.  Initial  Step 


As  discussed  previously,  these  are  complex  systems.  However,  at 
preb.nt,  the  available  data  seems  consistent  with  the  idea  that  the  most 
likely  mechanism  for  the  initial  decomposition  in  the  gaseous  phase  is  a 
combination  of  HOKO  elimination  at  lower  temperatures  and  N-NO2  cleavage  at 
higher  temp jratures,  with  the  crossover  point  occurring  at  about  300°C*”  or 
600°K.  This  is  based  on  thermochemical  estimates  of  temperature-variation 
of  relative  rates  for  the  possible  decomposition  steps, ^  on  comparison  of  the 
preliminary  best  values  (Table)  for  activation  energies  and  frequency  factors 
with  values  estimated27  for  the  individual  processes,  and  on  the  results  of 
very  low  *  essure  pyrolysis  studies1**  on  dimethylnltramine  decomposition. 


For  decomposition  in  the  pure  liquid  phase,  the  above  values  for  the 
activation  energies  for  decomposition  of  RDX  (47.8  kcal/mole)  and  HMX  (53-57 
kcal/mole)  seem  reasonably  close  to  the  values  estimated27  for  gas-phase  N-NO2 
cleavage  (46.2  kcal/mole),  although  the  value  for  HMX  does  seem  uncomfortably 
close  to  the  60  kcal/mole  estimated27  for  gas-phase  unassisted  C-N  cleavage 
and  is  also  associated  with  a  log  A  (sec-1)  value  of  1020  -  1021,  which  seems 
unexpectedly  high.  However,  without  more  detailed  knowledge  of  the  molecular 
and  electronic  structure  of  the  transition  state,  it  is  difficult  to  evaluate 
the  effect  of  such  factors  as  (a)  the  electrostatic  effect  of  the  other  N-NO2 
groupings  in  the  starting  HMX  or  FDX,  or  (b)  "solvent”  as  opposed  to  the  gas- 
phase)  effects  of  neat  liquid  HMX  and  RDX  on  the  observed  activation  energies. 

The  activation  energies  for  decomposition  of  HMX  and  RDX  in  solution  are 
difficult  to  interpret.  They  tend  to  be  around  40  kcal/mole,  which  is  several 
kcal/mole  lower  than  those  of  the  neat  liquids,  although  not  quite  as  low  as 
the  vapor-phase  values  (ca  30-35  kcal/mole).  The  reasons  for  this  are 
uncertain,  in  view  of  the  difficulty  in  evaluating  the  solvent  effects  on 


L.  Korsunskii,  F.  I.  Dubovitskil  and  G.  V.  Sitonina,  "Kinetics  of  the 
Thermal  Decomposition  of  N,N-Dimethylnitroamine  in  the  Presence  of 
Formaldehyde,"  Dokl.  Akad.  Nauk.  SSSR,  Vol.  174,  P.  1126,  1967 
(Engl.Transl. ,  p.  436). 


these  activation  energies.  One  possible  explanation  might  be  that  in  the  gas- 
phase  and  in  solvents  studied  to  date  the  first  step  is  HONO  elimination, 
while  in  the  neat  liquid  it  is  N-NO2  cleavage,  although  without  more 

information  it  is  difficult  to  suggest  why  the  "solvent  effect"  exerted  by 

neat  liquid  HMX  or  RDX  should  be  so  different  from  the  other  solvents.  Note, 
however,  that  the  studies  have  apparently  been  carried  out  in  a  temperature 
range  where  the  intrinsic  (gas  phase)  reactivities  toward  HONO  elimination  and 
N-NO2  cleavage  are  near  the  crossover  point;  thus  small  variations  in  solvent 

effects  could  conceivably  have  more  dramatic  effects  on  the  nature  of  the 

principal  pathway. 

Possibly  the  main  conclusion  to  be  drawn  from  the  activation  energies  is 
that  the  observed  values  for  HMX  and  RDX  decomposition  in  the  gaseous  and 
liquid  phases  seem  to  be  in  a  region  (ca  30-50  kcal/mole)  that  is  more 
characteristic  of  HONO  elimination  and  N-NO2  cleavage  than  of  the  higher 
values  estimated^'  for  C-N  cleavage,  although  it  is  difficult  to  rule  out  all 
possible  variations  on  the  general  theme  of  C-N  cleavage,  such  as  C-N  cleavage 
involving  two  or  more  bonds  simultaneously,  expecially  since  unexpectedly  low 
activation  energies  for  depolymerization  of  trioxane  to  formaldehyde  (a 
reasonable  model  for  depolymerization  of  RDX  to  N-nitroforminine)  ha  :  been 
interpreted  in  terms  of  a  concerted  ring  cleavage. However,  some  degree 
of  N-N  cleavage  must  be  capable  of  taking  place,  since  isotope-labeling 
studies  on  mixtures  of  un-  and  fully  N-15  labelled  HMX3^  show  some  s  ambling 
of  the  labelled  and  unlabelled  nitrogens.  The  scrambling  most  likely  takes 
place  during  the  primary  decomposition  step,  since  HMX  labelled  in  the  nitro 
group  with  N-15  decomposes  to  give  unscrambled  ^0  with  the  labelled  nitrogen 
still  attached  to  the  oxygen. 32,61  -j^e  mosC  straightforward  explanation  for 
the  scrambling3^  is  a  preliminary  equilibrium  involving  dissociation  and 
recombination  of  the  N-NO2  bond. 

However,  it  should  be  remembered  that  the  activation  energy  for  the 
decomposition  of  neat;  liquid  HMX  and  RDX  may  be  influenced  by  factors  other 
than  the  occurrence  of  a  single,  homogeneous  process.  For  example,  it  has 


59S.  W.  Benson  and  H.  E.  O'Neal,  "Kinetic  Data  on  Gas-Phase  Unimolecular 
Reactions,"  NSRDS-NBS  21,  p.  314ff.  For  sale  by  the  Superintendent  of 
Documents,  US  Printing  Office,  Washington,  DC  20402,  February  1970. 
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62  63 

been  suggested  ’  that  the  higher  values  observed  for  decomposition  of 
nitrate  esters  in  the  liquid  than  in  the  vapor  phase  may  be  due  to  a  decrease 
with  increasing  temperature,  in  the  inhibition  of  the  reaction  due  to 
recombination  of  NO^  and  alkoxy  radical  forged  in  the  primary  O-NO^  bond 
breaking  step.  This  decrease  was  suggested0  *  to  be  due  to  increased 
mobility  and  to  increased  tendency  for  the  to  undergo  further  reaction 
rather  than  recombination,  as  temperature  increases.  Like  the  nitrate  esters, 
the  nitramines  exhibit  higher  activation  energies  for  decomposition  in  the 
liquid  than  in  the  vapor  phase,  although  in  the  case  of  the  nitramines,  the 
gas-phase  activation  energies  and  frequency  factors  seem  more  characteristic 
of  HONO  elimination  than  of  N-NO2  cleavage.  However  the  above  "cage"  effect 
could  conceivably  be  raising  the  neat-liquid  activation  energies  for  HMX  and 
RDX  decomposition.  In  this  connection  it  seems  worth  mentioning  that 
nitroalkanes  appear  to  undergo  HONO  elimination  rather  than  C-Nitro  cleavage 
in  the  gas  phase,  at  least  up  to  temperatures  of  about  450°C. 

The  frequency  factors  are  somewhat  more  difficult  to  interpret, 
especially  the  condensed  phase  values.  The  frequency  factors  for  the  neat- 
liquid  decompositions  (Table)  PDX,  log  A  (sec-*)  *  18.7;  HMX  log  A(sec-1)  * 
19-22  seem  unexpectedly  high  for  a  first  order  reaction.  Although  this  is 
consistent  with  the  occurrence  of  bimolecular  follow-up  reactions  of  the  type 
discussed  previously,  it  should  be  remembered  that  little  if  anything  is 
known  about  energy  transfer  processes  in  liquid  HMX  and  RDX.  These  high 
frequency  factors  seem  inconsistent  with  bimolecular  initial  steps,  or  with 
anchimerically  assisted  C-N  cleavage,  since  ouch  mechanisms  would  be  expected 
to  lead  to  low  frequency  factors  due  to  the  high  degree  of  ordering  in  the 
transition  states. 

C.  Unimolecular  Follow-Up  Steps 

Some  possible  unimolecular  follow-up  steps  have  been  discussed 
previously;  these  include  (a)  ring  opening  by  0-cleavage  reactions,  and 
elimination  of  N-Nitroforroimine  (H2ONNO2)  from  the  ring-opened  open-chain 
intermediates  by  further  3-cleavage  reactions;  and  (b)  further  reaction  of 
H2ONNO2,  most  likely  to  give  formaldehyde  and  ^0  (both  of  which  are  known 
(see  for  example  References  1,  7-9,  43,  44  and  51)  to  be  formed  in  large 
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amounts  In  HMX  and  RDX  decomposition) ,  although  at  higher  temperatures  it  is 
possible  that  N-N  cleavage  of  N-Nitroformimine  *"o  NO2  and  I^CN  might  also 
occur. 

D.  Bimolecular  Follow-up  Steps 

As  discussed  earlier, ^  possible  bimolecular  follow-up  steps  include  (a) 
abstraction  of  hydrogen  from  starting  HMX  and  RDX  by  free  radicals  formed  in 
the  reaction,  and  (b)  attack  on  nitro  oxygen  by  free  radicals  formed  in  the 
reaction,  possibly  followed  by  N-0  or  N-N  cleavage  of  the  resulting 
oxynitroxlde  radical  to  give  a  nitrosoamine  or  a  denitro  nitrogen-centered 
radical.  Available  literature  analogies  and  thermo-chemical  estimates^  are 
consistent  with  the  idea  that  these  are  both  reasonable  follow-up  steps. 
Possibly  steps  such  as  these,  involving  radicals  formed  either  (a)  directly 
from  HMX  or  RDX  decomposition  or  (b)  by  reaction  between  primary  products  such 
as  Nf>2  or  ^0=0,  are  at  least  partly  responsible  for  the  high  frequency 
factors  and  apparent  autocatalysis  of  nitramine  decomposition  referred  to 
above. 


More  work  is  needed  on  all  aspects  of  HhK  and  RDX  decomposition 
mechanisms.  While  the  above  and  similar  mechanisms  seem  reasonable  and 
onsistent  to  the  present  writer,  the  individual  mechanisms  are  not  firmly 
stablished  and  in  fact  are  at  present  test  considered  only  as  a  basis  for 
further  discussion.  The  formation  of  other  known  products  of  HMX  and  RDX 
decomposition,  such  as  N7  and  NO,  can  be  explained  by  adaptations  or 
combinations  of  the  above  mechanisms,  products  such  as  00,  CD^t  anc*  ^0, 
could  be  formed  by  reaction  among  primary  products,  for  example  NO2  and  H?C0. 

E.  Extension  to  Combustion  Conditions 


Since  combustion  takes  place  at  much  higher  temperatures  and  pressures 
than  encountered  in  the  low-temperature  decomposition  studies  summarized  here, 
it  seems  appropriate  to  discuss  the  possible  effect  of  high  temperature, 
pressure  and  heating  rate.  These  have  been  discussed  in  our  previous 
presentation. ^  Some  further  discussion  follows. 

When  the  decomposition  takes  place  at  increased  pressure,  there  might  be 
an  increase  in  the  importance  of  bimolecular  processes  relative  to 
unimolecular  processes;  this  might  have  an  important  effect  on  combustion 
behavior,  as  discussed  previously. ^  Another  possible  effect  of  pressure  might 
be  that  at  sufficiently  high  pressures,  the  vapor  phase  might  become  so 
compressed  that  the  environment  "seen"  by  a  single  molecule  might  resemble  a 
liquid  more  than  a  low-pressure  vapor  or  gas;  if  this  were  the  case  the 
decomposition  mechanism  might  well  resemble  the  liquid  decomposition  rather 
than  a  low-pressure  gas;  for  modeling  purposes,  possibly,  assumed  activation 
energies  and  frequency  factors  should  be  modified  accordingly. 

The  effect  of  temperature  on  reaction  mechanisms  and  products  could  also 
be  important;  as  temperature  rises,  reactions  with  high  activation  energy 
would  be  expected  to  accelerate  relative  to  those  with  lower  activation 
energies,  since  the  activation  energy  is  simply  the  slope  of  a  plot  of  rate 
constant  V3  reciprocal  temperature.  Two  examples  of  this  have  appeared 
earlier  in  the  present  report;  these  are  the  apparent  temperature  dependence 
of  the  relative  rate  of  N-NO2  cleavage  and  of  H0N0  elimination  in  the  gas- 
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ohase  decomposition,^  and  the  possible  tendency  for  cleavage  of  to 

K7CN  and  N09  to  oecome  more  important  relative  to  formation  of  1^0  and  H2CO, 
ai  temperature  rises.  Such  changes  in  relative  importance  of  various  chemical 
mechanisms  could  give  rise  to  temperature-and  heating-rate-dependent  changes 
in  decomposition  product  distribution  and  chemistry;  as  with  the  pressure 
effects,  these  might  well  be  very  important  in  modeling  and  understanding 
combustion  and  explosive  behavior. 

Another  possible  effect  of  high  heating  rate  or  high  temperature  might  be 
to  cause  the  first  step  of  the  decomposition  and  unimolecular  follow-up  steps 
to  become  faster  relative  to  bimolecular  follow-up  steps  because  of  the  higher 
activation  energy  characteristic  of  first-order  reactions,  and  since  more 
unimolecular  decomposition  might  be  expected  to  take  place  immediately,  with 
corresponding  decrease  in  opportunity  for  bimolecular  follow-up  steps 
involving  starting  HMX  or  RPX  molecules  or  early  intermediates.  This  effect 
might  also  cause  important  changes  in  product  distribution  or  chemistry  with 
increasing  temperature  or  heating  rate.  The  temperature  and  pressure  effects 
on  the  unimolecular/bimolecular  ratio  apparently  work  in  opposite  direction 
but  complete  compensation  is  of  course  not  assured. 

Tn  view  of  the  preceding  paragraphs,  possibly  the  most  useful  aspect  of 
thermal  decomposition  studies  at  low  temperatures  and  pressures  is  not  to 
provide  product  distributions  and  kinetic  parameters  that  can  be  applied 
directly  to  combustion  conditions,  hut  rather  to  elucidate  the  types  of 
chemical  decomposition  processes  involved,  includitg  minor  (at  low 
temperatures  and  pressures)  pathways  in  addition  to  the  principal  ones. 
Informed  extension  of  this  body  of  knowledge  to  combustion  conditions  could 
then  provide  the  basis  for  improved  understanding  and  control  of  combustion 
processes  and  operational  properties  such  as  stability,  sensitivity,  and 
burning  rate  behavior. 


VI.  WORK  NEEDED 

More  work  is  needed  in  a  number  of  areas.  Several  of  these  seem 
particularly  worth  mentioning,  among  them  the  following: 

1.  The  first  step  of  the  decomposition  of  PMX  and  RDX  is  still  not  well 
understood.  In  addition  to  an  improved  understanding  of  the  nature  of  the 
first  step  itself  (for  example,  is  it  HOMO  elimination,  N-NOj  cleavage,  C-N 
cleavage,  or  something  else?)  there  is  a  need  for  a  better  understanding  of 
the  effect  of  temperature,  pressure  and  state  of  aggregation  (solid,  liquid, 
vapor,  as  well  as  solvent  effe^jjs^  Time-resolved  laser  spectroscopic 
techniques  may  be  helpful  here  ’  at  least  for  the  vapor-phase 
decomposition. 


66B.  H.  Rockney  and  E.  R.  Grant,  "Resonant  Multiphoton  Ionization  Detection  of 
the  NC>2  Fragment  from  Infrared  Multiphoton  Dissociation  of  CH2NP2."  Chem. 
Phys.  Lett. ,  Vol.  7S,  p.  15,  1981. 

^K.  E.  Lewis,  D.  F.  McMillen  and  D.  M.  Golden,  "Laser-Powered  Homogeneous 
Pyrolysis  of  Aromatic  Nitro  Compounds,"  J.  Phys.  Chem.,  Vol.  84,  p.  226, 
1980. 
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2.  There  seems  to  be  very  little  work  available  on  the  gas-phase 
decomposition  of  HMX  and  RDX.  There  is  a  need  for  more  understanding  of  this 
with  regard  to  both  the  kinetics  and  the  product  distributions,  especially 
since  the  gas-phase  decomposition  would  most  reflect  the  intrinsic  behavior  c 
the  molecules  themselves;  behavior  in  the  liquid  and  solid  phases  could  then 
be  regarded  as  variations  on  the  vapor-phase  behavior,  induced  by  the 
environment  set  up  by  neighboring  molecules.  In  addition  to  techniques  such 
as  those  described  in  References  16  and  66,  laser-powered  homogeneous 
pyrolysis**'  might  well  be  capable  of  providing  information  on  this  point. 

Note  that  care  should  be  taken  to  work  under  conditions  such  that  the  results 
are  not  influenced  by  wall  reactions. 

3.  There  is  a  need  for  improved  understanding  of  the  solid-ph.-se 
decomposition.  In  view  of  the  discussions  of  Cosgrove  and  Owen^~^'t  and  of 
Batten,  '  it  seems  clear  that  in  studying  the  behavior  Of  HMX  and  RDX  below 
their  melting  points,  care  should  be  taken  to  separate  the  true  solid-state 
decomposition  from  concurrent  sublimation  followed  by  gas-phase  decomposltioi 
and  to  understand  the  relative  roles  of  gas,  liquid,  and  solid-state 
decomposition;  see  the  discussion  on  pages  31-35  of  the  present  report. 

4.  There  is  also  a  need  for  Improved  understanding  of  the  effect  of 
temperature  and  pressure  on  the  decomposition  of  HMX  and  RDX.  In  extending 
low-temperature,  low-pressure  results  to  combustion  conditions,  it  should  be 
remembered  that  the  Arrhenius  parameters,  which  are  usually  used  to  carry  oul 
these  extensions,  may  themselves  depend  on  temperature  and  pressure;  improve* 
understanding  of  these  variations  would  enhance  the  reliability  of  the 
extension. 

5.  There  is  a  need  for  Improved  understanding  of  the  nature  and 
occurrence  of  autocatalysis  and  autolnhlbltlon  in  the  thermal  decomposition  c 
HMX  and  RDX.  This  is  especially  true  in  view  of  apparent  discrepancies 
between  literature  statements  concerning  acceleration  or  inhibition  of  HMX  at 
RDX  decomposition  by  known  products  such  as  formaldehyde  and  oxides  of 
nitrogen.^ .46,49,68  gee  ai80  the  discussion  above  under  ’’Autocatalysis"  an< 
references  cited  therein. 


6.  There  is  also  a  need  for  more  research  Into  the  effect  of  pressure  i 
HMX  and  RDX  decomposition,  with  regard  to  both  kinetics  and  product 
distributions.  One  study,  ’  carried  out  at  pressures  up  to  1000  psi  (0.067 
kbar),  resulted  in  activation  energies  for  HMX  decomposition  which  decreas  id 
from  55.9  kcal/mole  at  atmospheric  pressure  to  41  8  kcal/mole  at  1000  psi 
(0.067  kbar)  of  argon,  although  the  activation  energy  for  RDX  decomposition 


D.  F.  Dehenham  and  A.  J.  Owen,  "The  Thermal  Decomposition  of  1,3,5- 
Trinitrohexahydro-1 ,3 ,5-triazine  (RDX)  in  1 ,3,5-Trinitrohenzene,"  Symp. 
Chem.  Probl.  Connected  Stab.  Explos.  (Proc.),  Vol.  4,  p.  201,  1976.,  Pub. 
1977.  \ 


*D.  A.  Flanigan  and  B.  B.  Stokes,  "HMX  Deflagration  and  Flame 
Characterization.  Volume  I.  Phase  II  Ni^ramine  Decomposition  and 
Deflagration  Characterization,"  Thiokol  Corporation,  Huntsville  Division, 
Huntsville,  AL,  Report  AFRPL-TR-79-94,  October  1980  (AD-B053  058L). 
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was  reported  to  be  independent  of  pressure  from  atmospheric  pressure  to  500 
psi  (0.03  kbar) .  In  a  study  of  time-to-explosion  versus  temperature  for  a 
number  of  explosives  including  HMX,^  it  was  found  that  "activation  energy" 
appeared  constant  from  <1  kbar  to  50  kbar;  it  was  suggested  that  the 
decomposition  mechanism  was  independent  of  pressure.  There  has  also  been  a 
report^  that  the  activation  energy  for  decomposition  of  RDX  was  ca  15 
kcal/mole  lower  when  measured  by  closed-pan  than  by  open-pan  DSC.  Pressure 
effects  on  product  distributions  from  HMX  and  RDX  decomposition  have  been 
summarized  briefly.^  In  view  of  the  differences  between  experimental 
conditions  and  results  reported  in  these  studies , ^ » 69—7 1  c^ear  t^e 

pressure  effects  on  HMX  and  RDX  decomposition  are  far  from  well  understood, 
and  much  more  work  in  this  area  is  required. 

VII.  MOTHS  ADDF.D  IN  PROOF 

A.  Note  Added  in  Proof  No.  1  While  this  report  was  being  proofread,  the 
author  became  aware  of  a  review  article^  in  the  Russian  literature;  although 
as  this  is  written  this  review  probably  will  not  be  available  in  English 
translation  for  several  months,  it  seems  worth  while  to  call  it  to  the 
attention  of  anyone  reading  the  present  report. 

B.  Note  Added  in  Proof  No.  2  Brill  and  Karpowicz^ have  recently 
suggested  that  the  decomposition  reactions  of  condensed-phase  HMX  and  RDX  may 
be  governed  by  release  of  intermolecular  interactions  occurring  prior  to 
unimolecular  decomposition,  rather  than  by  the  decomposition  steps 
themselves.  This  is  based  on  large  (ca  40-50  kcal/mole)  activation  energies 
observed  for  phase  transitions  among  the  solid  polymorphic  forms  of  HMX;  Brill 
and  Karpowicz  feel  that  the  large  act ivat ion  energies  for  these  phase 
transitions  (much  larger  than  the  net  energies  of  reaction  for  the  same 
transitions)  represent  energy  needed  for  "freeing-up"  of  HMX  or  RDX  molecules 


E.  R.  Lee,  R.  H.  Sanborn  and  H.  D.  Stromberg,  "Thermal  Decomposition  of  High 
Explosives  at  Static  Pressures  10-50  Kilobars,"  Proceedings  5th  Symposium 
(International)  on  Detonation,  1970,  p.331-337,  for  sale  by  the  Superintendent 
of  Documents,  U.S.  Government  Printing  Office,  Washington,  Du  20402. 

*M.  A.  Schroeder,  "Critical  Analysis  of  Nitramine  Decomposition  Data: 

Product  Distributions  from  HMX  and  RDX  Decomposition,"  Proceedings  18th 
JANNAF  Combustion  Meeting,  Pasadena,  CA,  Vol.  2,  p.  395,  CPIA  Publication 
347,  October  19-23,  1981. 

^F.  I.  Dubovitskii  and  B.  L.  Korsunskii,  "Thermal  Decomposition  Kinetics  of 
N-Nitro  Compounds,"  Usp.  Khim.,  Vol.  50,  p.  1828,  1981,  cited  ir.  Current 
Contents,  Physical,  Chemical  and  Earth  Sciences,  Vol.  22,  p.144.  No.  2, 
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T.  B.  Brill  and  R.  J.  Karpowicz,  "Solid  Phase  Transition  Kinetics  The  Role 
of  Intermolecular  Forces  in  the  Condensed- Phase  Decomposition  oe  Octahydro- 
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Transition  Kinetics  to  the  Properties  of  HMX,"  AIAA  Journal,  Vol.  21, 
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42 


which  presumably  must  take  place  before  decomposition.  The  magnitudes  of  the 
energies  involved  are  such  as  to  raise  at  least  the  possibility  that  the 
freeing-up  of  the  molecules  may  be  the  rate-determining  step,  at  least  for  the 
solid-state  decomposition. 

Note  that  this  treatment  depends  critically  on  the  assumption  that  the 
rate  constants  measured  in  References  73  and  74  and  attributed  to  solid-solid 
phase  transitions  in  HMX,  are  actually  rate  constants  for  true  first-order 
processes  and  hence  that  the  activation  energies  are  also  true  activation 
energies  for  a  specific  first-order  process,  not  merely  effective  temperature 
coefficients  of  reaction  (see  Reference  24)  for  a  complicated  assembly  of 
processes. 

In  addition,  Fifer7^  has  published  a  more  recent  and  accessible 
discussion  than  reference^  Qf  the  cage  effect  hypothesis.  The  basic  Idea  of 
this  hypothesis  is  that  the  relative  ordering  of  the  rates  and  Arrhenius 
parameters  for  the  gas,  liquid  and  solid-state  thermal  decomposition  of 
nitramines  and  nitrate  esters  (Gas  phase  rates  faster  than  liquid,  which  in 
turn  are  faster  than  solid-state  rates  of  decomposition;  Arrhenius  parameters 
follow  the  opposite  ordering)  is  explainable  in  terms  of  greater  confinement 
of  radicals  (such  as  NU2  and  denitro-HMX/RDX)  formed  in  the  initial  step  of 
decomposition.  This  is  in  agreement  with  recent  work‘d  on  the  thermolysis  of 
l  ,2-diphenylethane. 

Note  that  this  hypothesis  applies  only  to  radical-producing  reactions  and 
not,  for  example,  to  HONO  elimination.  However,  even  if  the  gas-phase 
reaction  were  found  ultimately  to  proceed  via  some  other  mechanism,  for 
example  HONO  elimination,  the  cage-effect  hypothesis  might  well  still  be 
useful  to  understand  quantitative  differences  between  the  observed  liquid- 
phase  decomposition  kinetics  and  parameters,  and  values  reliably  estimated  for 
the  hypothetical  gas-phase  radical  pathways. 

Furthermore  it  was  argued7^  that  a  practical  implication  of  the  cage 
effect  hypothesis  was  that  at  higher  temperatures  the  condensed  phase  rates 
would  approach  those  for  the  gas  phase;  so  that  at  the  higher  temperature 
characteristic  of  combustion  it  might  be  a  reasonable  approximation  to  use  the 
gas  phase  kinetic  parameters  to  represent  the  condensed  phase  decomposition. 
The  present  writer  suggests  it  may  also  be  worthwhile  to  consider  the 
possibility  that,  at  the  higher  pressures  characteristic  of  combustion,  the 
density  of  the  gas  phase  might  be  increased  to  the  point  where  it  might 
resemble  the  liquid  more  than  what  would  be  thought  of  at  more  moderate 
pressures  as  a  vapor:  if  this  is  true,  the  system  might  exhibit  behavior  such 
that  it  would  be  the  vapor  phase  that  would  assume  Arrhenius  parameters 
characteristic  of  the  liquid,  rather  than  the  other  way  around. 


7^R.  A.  Fifer,  "Cage  Effects  in  the  Thermal  Decomposition  of  Nitramines  and 
Other  Energetic  Materials,”  Proceedings  19th  JANNAF  Combustion  heeting, 
Greenbelt,  MD;  Vol.  I,  CPIA  Publication  366,  p.  311-319,  October  1982. 

76S.  E.  Stein,  D.  A.  Robaugh,  A.  D.  Alfieri  and  R.  E.  Miller,  "Bond  Homolysis 
in  High-Temperature  Fluids,"  J.  Amer.  Chem.  Soc..  Vol.  104.  p.  6567-70. 
1982. 
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At  the  very  least,  both  of  the  effects^-' ^  described  in  the  above 
paragraphs  will  have  to  be  batter  understood  than  at  present  before  Arrhenius 
parameters  applicable  to  combustion  conditions  can  be  reliable  derived  from 
studies  of  decomposition  under  more  moderate  temperatures  and  pressures.  It 
will  also  of  course  be  necessary  to  understand  the  role  of  such  factors  as 
dielectric  constant  of  the  medium,  and  other  types  of  solvent  effects  and 
lntermolecular  interactions, 

C.  Note  Added  in  Proof  No.  3  It  also., seems  worthwhile  to  mention  an 
interesting  paper  in  which  Burov  and  Nazin^  have  very  recently  reported  some 
kinetic  measurements  on  gas  phase  thermal  decomposition  of  several  nitramlnes 
Including  HMX  and  RDX.  Other  compounds  studied  included  N,N-dinitropiperazine 
and  1 ,3-dinitrazacyclopentadiene.  The  decompositions  of  all  of  these 
compounds  except  HMX  were  studied  In  the  presence  of  NO,  addition  of  which 
noticeably  affected  the  kinetics,  resulting  in  lower  rates  by  factors  of  2  to 
4  at  200°C,  as  well  as  noticeable  changes  in  Arrhenius  parameters;  thus  chain 
reactions  or  autocatalysis  are  apparently  of  some  importance  here.  The 
decompositions  were  studied  manometrlcally,  at  pressures  of  0.05-0.1  mm  Hg. 

Gas  phase  and  solution  kinetic  measurements  on  the  above  nitramlnes  as  well  as 
dimethyl  and  diethyl  nitramlnes  and  N-nitroplperidlne  were  summarized  in  a 
table. " 

The  authors^  make  the  interesting  suggestion  that  differences  in 
decomposition  rates  of  these  nitramlnes  are  determined  by  the  configuration  of 
the  amino  nitrogens,  the  more-rapidly  decomposing,  low-activation-energy 
compounds,  including  RDX,  being  those  with  tetrahedral  or  pyrlmidal  nitrogens, 
while,  the  slower-decomposing,  higher  activation  energy  compounds  have  trigonal 
coplanar  nitrogens.  This  will  be  discussed  further  in  a  final,  wrap-up 
portion  of  the  present  review,  which  will  be  devoted  to  chemical  and  physical 
mechanisms  of  decomposition. 

It  also  seems  worth  mentioning  that  the  authors  of  Reference  77  are  of 
the  opinion  that  the  N-NtU  bond  dissociation  energy  of  nitromethane  is  about 
40  kcai/mole,  giving  as  their  authority  a  Russian  book'8  which,  as  far  as  the 
present  writer  is  aware,  is  not  available  in  English  translation.  This  value 
is  considerably  lower  than  the  46.2  kcal/mole  estimated  by  Shaw  and  Walker  for 
dlmethylnitramine  and  applied  also  to  HMX.  Korsunskii  and  Dubovltskii'^  also 
seem  to  feel  that  the  Jf-N02  bond  dissociation  energy  for  simple  nitramlnes  is 
around  40  kcal/mole.  This  discrepancy  should  be  checked  into  further,  in  view 
of  the  obvious  importance  of  an  accurate  knowledge  of  the  K-N02  dissociation 
energies  of  dlmethylnitramine  and  other  nitramlnes  to  the  decomposition 
chemistry  of  these  compounds.  In  view  of  the  possible  Importance  of  wall 
reactions  and  autocatalysis,  and  of  decomposition  pathways  such  as  HONO 


^Yu.  M.  Burov  and  G.  M.  Nazin,  "Influence  of  Structure  on  Rate  of 
Decomposition  of  Secondary  Nitramlnes  in  the  Gas  Phase,"  Klnetika  i  Katallz, 
Vol.  23,  No.  L,  pp.  12-17,  1982  (English  Translation,  p.  5-10). 
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V.  A.  Dedvedev,  V.  K.  Potapov  and  Yu.  S.  Khodeev,  "Chemical  Bond  Rupture 
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elimination,  this  knowledge  should  be  obtained  from  sources,  such  as 
thermochemical  studies,  that  are  independent  of  kinetic  measurements  on  the 
decomposition  reactions  of  the  nitramines. 

It  seems  worth  noting  that  Reference  77  says  that  Shaw  and  Walker's 
estimated  parameters  for  N-NO2  bond  dissociation  were  taken  from  a  private 
communication  concerning  DMNA.  The  value  of  E  was  stated  to  have  been  adopted 
without  change,  and  log  A  to  have  been  multiplied  by  4.  However,  an 
examination  of  Shaw  and  Walker's  paper  (Reference  27  of  the  present  report) 
shows  that  they  are  based  on  heats  of  formation  and  frequency  factors  given  in 
References  12-14  and  21  of  Reference  27  of  the  present  report;  these 
’•eferences  are  not  to  any  private  communications  but  to  two  journal  articles 
and  to  two  standard  compilations  of  data.  In  particular  the  N- NO 2 
dissociation  energy  is  based  on,  in  addition  to  (013)2  NNO2,  heats  of 
formation  for  (CHj^N*  and  NfU,  which  are  the  same  as  those  given  in  standard 
compilations  of  data.®®*®*  Shaw  and  Walker  also  cite  a  number  of  reported 
values  for  the  N-NO2  bond  dissociation  energy,  nearly  all  of  which  are  higher 
than  40  kcal/mole,  ranging  from  47  to  66  kcal/mole.  Thus  the  46.2  kcal/mole 
estimated  by  Shaw  and  Walker  for  the  N-NO2  dissociation  energy  of 
dimethylnitramine  seems  to  be  reasonably  well-founded  in  the  literature. 
However  it  would  be  most  helpful  to  have  access  to  the  reference7”  cited77  as 
supporting  an  upper  limit  of  170  kJ/mole  (40.6  kcal/mole)  for  the  N-NO2  bond 
energy  of  dimethylnitramine;  expecially  since  additional  measurements  might 
lead  to  revision  of  the  heats  of  formation  on  which  the  figure  of  46.2 
kcal/mole  was  based.  The  frequency  factor  (log  A  ■»  15.8)  suggested  by  Shaw 
and  Walker27  was  based  on  an  average  of  values  from  decomposition  of 
nltromethane  and  ^0^;  possibly  more  elaborate  thermochemical  calculations 
would  be  helpful  in  arriving  at  a  firmer  value. 

Activation  energies  and  frequency  factors  (log  A(sec“*))  reported  in 
Reference  77  were  as  follows:  for  RDX,  146.6  kJ/mole  (35.04)  kcal/mole)  and 
13.5  without  NO;  167.6  kj/mole  (40.06  kcal/mole)  (see  footnote,  p.  6  of 
Reference  77)  and  15.6  with  NO.  For  HMX,  the  values  were  165.5  kj/mole 
(39.56)  kcal/mole)  and  14.2  without  added  NO.  N,N’-dinitropiperazine  gave  160 
kj/mole  (38.2  kcal/mole)  and  13.6  alone,  but  151.2  kj/mole  (36.14  kcal/mole) 
and  12.0  in  the  presence  of  NO,  while  1 ,3-dinitrazacyclopentane  gave  146.6 
kj/mole  (35.04  kcal/mole)  and  13.5  alone  and  169.3  kj/mole  (40.46)  kcal/mole) 
and  15.6  in  the  presence  of  NO. 

When  the  results  for  RDX  are  plotted  alongside  those  given  in  Figure  2, 
they  do  not  fall  on  the  same  line  as  the  other  results  plotted  in  Figure  2, 
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but  lie  on  a  parallel  line,  about  0.5  log  unit  below  those  plotted  in  the 
figure.  However,  they  no  not  seem  to  fit  the  other  points  on  the  graph  as 
well  as  the  data  of  Belyaeva  et  al.,  judging  from  the  degree  of  scatter  when 
the  line  defined  by  the  points  of  Reference  77  is  extended  through  the  other 
points. 


There  also  seems  to  be  a  series  of  typographical  errors  in  Table  1, 
column  6  of  Reference  77;  the  values  given  do  not  mate  »  those  calculated  from 
the  values  of  E  and  log  A  (sec-*)  in  columns  4  and  5.  The  present  writer 
recalculated  the  numbers  in  column  6  of  Table  1  of  reference^  and  obtained 
the  following  values:  line  1,  0.268  instead  of  4;  line  2,  1.67  instead  of 
0.158;  line  3,  5.0  instead  of  0.525;  line  4,  721  instead  of  0.676;  line  5, 

2.12  instead  of  0.2;  line  6,  31.2  instead  of  2.9;  line  7,  22.5  instead  of 
2.24;  line  8,  8.73  instead  of  0.87;  line  9,  530  instead  of  50;  line  10,  666 
instead  of  63;  line  11,  207  instead  of  20;  line  12,  81.3  instead  of  8;  line 
13,  1.7  instead  of  0.56;  line  14,  5.20  instead  of  0.56;  line  15,  255  instead 
of  2500;  line  16,  8.50  instead  of  0.8;  line  17,  0.631  instead  of  0.063;  line 
18,  8.64  instead  of  0.83;  line  19,  2.03  instead  of  0.22;  line  20,  24.5  instead 


of  2.5;  line  21,  207  instead  of  20;  line  22,  81.3  instead  of  8. 
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p.  59  (summarized  in  Chem.  Abstr.,  Vol.  90,  206758),  1978. 

D-28  A.  E.  Axworthy,  J.  E.  Flanagan,  and  J.  C.  Gray,  "Interaction  of  Reaction 
Kinetics  and  Nitramine  Combustion,"  Rocketdyne  Division,  Rockwell 
International,  Canoga  Park,  CA,  Report  AFATL-TR-80-58 ,  May  1980  (AD-B052 
861L) . 

D-29  M.  Boudart,  "Kinetics  of  Chemical  Processes,"  Prentice-Hall,  New  York, 
Chapter  6,  1968. 
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APPENDIX  A 


REPORTED  ARRHENIUS  PARAMETERS  FOR  THERMAL  DECOMPOSITION 
OF  HMX  AND  RDX  IN  THE  VAPOR  PHASE 

Available  values^- for  gas-phase  decomposition  of  HMX  and  RDX  are 
summarized  in  Table  A-l,  together  with  temperature  ranges,  kinetics,  etc.  The 
activation  energy  and  frequency  factor  for  gas  phase  decomposition  of  HMX  have 
been  determined  by  isothermal  DSCA-*,A  (see  discussion  accompanying  Appendix 
B,  p.  70-72),  by  very-low-pressure  pyrolysisA“2  and  by  a  methodA-3  involving 
combined  manometric  study  of  gas-  and  condensed-phase  decomposition  in  the 
very  early  stages  of  the  reaction.  The  DSC  valuesA~*  are  much  higher  than 
those  determined  by  the  other  methods. A-2 ,A-3  Both  DSC  and  VLPP  data  for 
HMXA-*’2  show  considerable  scatter,  but  the  VLPPA-2  values  seem  more  reliable 
due  to  the  larger  temperature  range  (135°  as  compared  to  12°);  thus,  the 
activation  energy  for  gas-phase  HMX  decomposition  in  the  temperature  range  ca_ 
250-400#C  would  appear  to  be _ca  33  kcal/mole,  with  a  frequency  factor  of  log  A 
(sec-*)  ca_  12.8  although  the  authors  point  out  that  their  data  are  too 
scattered  \o  completely  rule  out  that  45.50  kcal/mole  that  would  signify 
initial  N-NO2  bond  scission. A-*  This  value  of  £a^  33  kcal/mole  is  in  agreement 
with  the  manometric  values. A~ 3 

The  DSC  valuesA-*»A  for  RDX  decomposition  in  the  vapor  phase  were 
measured  over  a  wider  temperature  range  (51°)  than  for  HMX  (12°),  and  the  plot 
shows  considerably  less  scatter;  thus  the  RDX  values  measured  by  DSC  seem  much 
more  reliable  than  the  HMX  values  measured  by  DSC.  For  further  discussion  of 
DSC  studies  on  vapor-phase  HMX  and  RDX  decomposition,  see  Appendix  B. 

Since  HMX  and  RDX  might  be  expected  to  decompose  by  similar  mechanisms, 
it  is  interesting  to  note  the  general  agreement  between  the  manometricA-3  and 
VLPPA-2  values  for  HMX  and  the  values  for  RDX.A-*,3,A 


A_1R.  N.  Rogers  and  G.  W.  Daub,  "Scanning  Calorimetric  Determination  of  Vapor- 
Phase  Kinetics  Data,"  Anal.  Chem. ,  Vjl.  45,  p.  596,  1973. 

A-2D.  F.  McMillen,  J.  R.  Barker,  K.  E.  Lewis,  P.  L.  Trevor  and  D.  M.  Golden, 
"Mechanisms  of  Nitramine  Decomposition:  Very  Low-Pressure  Pyrolysis  of  HMX 
and  Diraethylnitramine,"  Final  Report  on  SRI  Project  PYU  5787,  Supported  by 
Department  of  Energy  through  Lawrence  Livermore  Laboratories,  Contract  No. 
EY-7 6-C-03-C 115,  18  June  1979. 

A_3M.  S.  Belyayeva,  G.  K.  Klimenko,  L.  T.  Babaytseva,  P.  N.  Stolyarov, 

"Factors  Determining  the  Thermal  Stability  of  Cyclic  Nitramines  in  the 
Crystalline  State,"  FTU-XD( RS)T- 1209-78  (AD-B033  572L). 

A  AR.  M.  Potter,  G.  J.  Knutson,  and  M.  F.  Zimmer,  Corapat.  Propellants  Explos. 
Pyrotechnics  Plast.  Addit.,  Conf.  Pub.  1975,  II-D.,  12  pp. ,  Chem.  Abstr., 
Vol.  86,  14241 lj  ,  1974. 
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As  discussed  In  the  main  body  of  the  report,  reactions  on  the  container 
surface  probably  are  not  influencing  these  values  too  much,  although  more  work 
on  this  point  would  help  since  wall  reactions  do  seem  to  influence  dimethyl- 
nitr amine  decomposition. 

Tables  A-2  and  A-3  show  the  actual  data  points  reported  by  the  authors  of 
References  A-l-4  these  were  analyzed  statistically  as  described  in  the  main 
report  (p.  19-26).  Recommended  values  are  given  in  the  table  in  the  main 
report,  and  the  points  are  plotted  in  Figures  2  and  4  of  the  main  report. 


rj)X  230-50  503-23  first  order  1I.70±0.78  29. 6i±1.85  12418  A-4  Isothermal  DSC 
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APPENDIX  B 


REPORTED  ARRHENIUS  PARAMETERS  FOR  DECOMPOSITION 
OF  NEAT  LIQUID  AND  DISSOLVED  HMX 

Available  values®^-*®  for  liquid-phase  decomposition  of  HMX  are 
summarized  in  Table  B-l,  and  solution-phase  values  in  Table  B-2,  together  with 
temperature  ranges,  kinetics,  etc. 

The  best  literature  set  of  Arrhenius  parameters  for  thermal  decomposition 
of  neat  liquid  HMX  (Table  Bl)  seems  to  be  the  52.7  kcal/mole  and  log  A(sec”*) 
of  19.7  reported  by  Robertson;  this  is  based  on  the  long  temperature  range 


®~^A.  J.  B.  Robertson,  "The  Thermal  Decomposition  of  Explosives.  Part  II 
Cyclotrimethylenetrinitramine  and  Cyclotetramethylenetetranitramine,” 

Trans.  Faraday  Soc«,  Vol.  45,  p.  85,  1949. 

®“^R.  h.  Rogers  and  E.  D.  Morris,  Jr.r  "On  Estimating  Activation  Energies  with 
a  Differential  Scanning  Calorimeter,"  Anal.  Chem.,  Vol.  38,  p.  412,  1966. 

®~^R.  N.  Rogers  and  L.  C.  Smith,  "Estimation  of  Preexponential  Factor  from 
Thermal  Decomposition  Curve  of  an  Unweighed  Sample,"  Anal.  Chem.,  Vol  39, 
p.  1024,  196# . 

®“^E.  K.  Rideal  and  A.  J.  B.  Robertson,  "The  Sensitiveness  of  Solid  High 
Explosives  to  Impact,"  Proc.  Roy.  Soc.,  Vol.  195A,  p.  135,  1948. 

®”^R.  N.  Rogers,  "Differential  Scanning  Calorimetric  Determination  of  Kinetic 
Constants  of  Systems  chat  Melt  with  Decomposition,"  Thermochimia  Acta, 

Vol.  3,  p.  437,  1972. 

®“^B.  B.  Goshgarlan,  “The  Thermal  Decomposition  of  Cyclotrimethylene¬ 
trinitramine  (RDX)  and  Cyclotetramethylenetetranltraralne  (HMX),"  AFRPL-TR- 
78-76  (AI>-B032  275L). 

®~^D.  A.  Flanigan  and  B.  B.  Stokes,  "HMX  Deflagration  and  Flame 

Characterization.  Volume  I.  Phase  II  Nitrsmine  Decomposition  and 
Deflagration  Characterization,"  Thiokol  Corporation,  Huntsville  Division, 
Huntsville,  AL,  Report  AFRPL-TR-79-94,  October  1980  (AD-B053  058L). 

®~®R.  J.  Powers,  "Hazards  Due  to  Precorab^istion  Behavior  of  High  Energy 

Propellants,”  AFATL,  Private  Communication,  June  1980,  of  data  presented  at 
JANNAF  Workshop,  Eglin  AFB,  April  1980. 

B  ft  l 

Yu.  Ya.  Maksimov,  Thermal  Decomposition  of  Hexagen  and  Octogen,"  Tr.  Mosk. 

Khim-Tekhnol ,  Inst.  No  53,  p.  73,  1967!  see  Chem.  Abstr.,  Vol.  68, 
p.  41742r.  Translated  by  H.  J.  Dahlby;  Los  Alamos  Report  LA-TR-68-30,  Los 
Alamos,  1968.  1 

®“l°Maksiraov  (presumably),  work  summarized!  in  K.  K.  Andreyev,  "Thermal 

Decomposition  and  Combustion  of  Explosives,"  FTD-HT-23-1329-68,  p.  87, 
p.  130  (AI>-693  600).  1 
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covered  (44°);  on  the  fact  that  the  absence  of  any  effect  of  sample  size 
(varied  from  l-7mg)  suggests  that  self  heating  is  not  important;  and  on  the 
low  degree  of  scatter  in  the  Arrhenius  plot.  However  some  of  the  other 
authors  do  not  state  temperature  range  or  provide  enough  information  to  enable 
the  degree  of  scatter  to  be  evaluated.  On  reading  Robertson's  paper  and 
experimental  reference  cited  therein,  no  obvious  sources  of  systematic  error 
were  evident  to  this  writer. 

Most  of  the  values  1  of  A  and  Ea,  measured  at  or  below  atmospheric 
pressure,  given  in  Table  B1  for  decomposition  of  neat  liquid  HMX  agree  within 
a  few  kcal/mole  with  those  of  Robertson.  The  exceptions  are  the  values  of 
References  B-6  and  B-8.  However,  these  points  were  measured  over  the  shortest 
temperature  ranges  in  Table  B-l  and  it  is  difficult  to  rule  out  the 
possibility  of  self  heating,  since  Reference  B-6  states  (p.10)  that  the  sample 
size  was  2  rag,  but  it  is  stated  on  page  34  that  ‘"sample  weight  greater  than 
1 .Omg  increased  the  isothermal  temperature  during  HMX  decomposition, 
apparently  due  to  excess  thermal  evolution  beyond  the  capabilities  of  the 
instrument  to  adequately  regulate  a  constant  temperature." 

Reference  B-8  used  sample  sizes  on  the  order  of  1.5  mg,  so  self-heating 
may  not  be  quite  as  bad.  Robertson  used  sample  sizes  of  1-7  mg  and  found  no 
eff  :ct  of  mass  ie.,  apparently  no  self-heating.  Possibly  this  is  because  in 
Lh«  2-cra  glass  bulb  or  copper  oven  apparatus  used  by  Robertson,  the  sample 
would  be  more  spr.  ad  out  than  in  a  millimeter  s  zed  DSC  pan. 

Judging  from  the  results  in  Table  B-l,  activation  energies  and  frequency 
factors  seem  to  be  independent  of  cype  and  pressure  of  inert  gases  present  up 
to  1  atmosphere;  however,  judging  from  the  data  of  Reference  B-7  it  would 
appear  that  on  going  from  1  to  70  atmospheres,  the  activation  energy  decreases 
noticeably,  from  ca_  56  kcal/mole  at  1  to  ££  42  kcal  at  70  atmospheres.  These 
results  are  least  squares  averages  and  appi  ir  independent  of  particle  size. 

In  addition  to  the  data  given  in  Tables  B-l  and  B-2,  R.  N.  Rogers®^** 
furnished  several  points  for  decomposition  of  neat  liquid  HMX,  measured  by 
isothermal  DSC  with  correction  for  the  concurrent  gas-phase  reaction. 
Powers,®*®  on  the  other  hand,  in  his  studies  of  the  decomposition  of  neat 
liquid  HMX  by  isothermal  DSC,  felt  that  the  gas-phase  decomposition  was 
sufficiently  unimportant  to  be  neglected;  this  is  in  agreement  with 
Goshgarian' observation  that  the  decay  portion  of  HMX  decomposition  curves 
appeared  smooth  throughout  the  decomposition.  However  for  spread  ^ample  a 
gas-phase  "hump"  could  be  seen.  Thus  the  gas-phase  decomposition  seems  less 
important  for  HMX  than  for  RDX.  In  any  case.  Powers'  data  were  obtained  in 
the  initial  phase  of  the  liquid  decomposition,  where  the  vapor  phase  should  be 
proportionately  least  important;  thus  they  were  included  along  with 
Robertsons®**  and  Rogers'  gas-phase-corrected  data,  in  the  best  three  sets  of 
data  on  liquid-phase  HMX  decomposition. 

The  actual  data  points  obtained  by  Robertson,®”*  Rogers®"**  and  Powers®"® 
are  given  in  Table  B-3.  These  points  were  analyzed  statistically  as  given  in 


B- 11 

iAR.  N.  Rogers,  Los  Alamos  Scientific  Laboratories, 
February  1980. 


Private  Communication, 
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the  main  report,  and  plotted  In  Figure  3  of  the  main  report.  Note,  however, 
that  while  the  picture  of  concurrent  liquid  and  vapor-phase  reactions  being 
the  main  contributors  to  the  shape  of  the  DSC  curves  still  seems  to  hold  for 
RDX,  it  has  recently  begun  to  appear®"®*'*  that  the  situation  for  HMX  may  be 
more  complicated,  possibly  as  a  result  of  complications  due  to  presence  of  th 
non-volatile  residue.  Until  this  situation  is  understood  better  than  it  is  a 
present,  it  may  well  be  a  good  idea  to  just  use  Robertson's  manometrlc  data  a 
the  best  data  for  decomposition  of  liquid  HMX.  As  shown  in  Table  B-3,  this 
gives  Efl  »  53.1  kcal/mole  and  log  A  (sec-1)  ■  19.950. 

There  are  only  a  few  values  for  decomposition  of  HMX  in  solution.  The 
value  of  50.3  kcal/mole  in  Convalex-10  was  measured  by  dynamic  DSC,  and  there 
is  no  information  on  the  temperature  range  covered.  The  data  of 
Maksimov®"’’'  were  done  over  a  fairly  large  / 54°C)  temperature  range,  and  th 
degree  of  scatter  in  the  Arrhenius  plot  given  does  not  seem  unreasonable. 
Thus  it  does  not  seem  unreasonable  to  conclude  from  the  data  in  Tables  B-2  an 
B-l,  that  the  activation  energy  and  frequency  factor  for  decomposition  of  HMX 
in  an  inert  solvent  are  somewhat  lower  than  those  for  neat  liquid  HMX, 
although  the  data  for  HMX  in  solution  were  aot  reanalyzed  statistically  by  th 
present  writer. 
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APPENDIX  C 


REPORTED  AR> HENIUS  PARAMETERS  FOR  DECOMPOSITION 
OF  N.  AT  LIQUID  AND  DISSOLVED  RDX 

Available  kinetic  parameters'”*”^*  for  decomposition  of  RDX  in  the  neat 
liquid  phase  and  in  solution  are  summarized  in  Tables  C-l  and  C-2,  respectively. 
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Secondary  Explosives  RDX  and  HMX,"  Defense  Science  Journal,  Vol.  18, 
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C-^*P.  W.  M.  Jacobs  and  A.  R.  T.  Kureishy,  J.  Chem.  Soc.,  p.  4718,  1964. 


Due  to  the  long  temperature  range,  low  degree  of  scatter  in  the  Arrhenius 
plots,  and  care  to  avoid  such  sources  of  error  as  self  heating,  the  best  pair 
of  Arrhenius  parameters  for  the  decomposition  of  liquid  RDX  is  probably  that 
given  by  Robertsonc~* (Table  C-l). 

However,  there  remains  the  problem  of  vaporization  followed  by  gas-phase 
reactions.  Robertson's  other  set  of  datac-4  are  probably  in  error  as  Wilbyc-5 
quotes  Reference  C-4  as  saying  that  the  values  may  be  in  error  due  to  poor 
equilibrium  and  to  gas-phase  decomposition*  The  values  of  Bawn,  Hays  and 
Pollardc-’'  and  of  Robertson  and  SinclairC-/  are  probably  in  errorc-5  due  to 
short  temperature  range  and  large  sample  sizes  leading  to  self  heating.  Along 
with  Robertson's  values, the  values  of  G.  K.  Adams''-  were  accepted  by 
Wilby  as  the  best  values.  However,  examination  of  Table  1  of  Adams' 
reportc-J  reveals  that  while  the  values  for  the  experimentally  determined 
logjgkex  and  the  calculated  log.Q  'cc  match  each  other  they  did  not  match 
values  calculated  from  the  stated  Arrhenius  parameters,  and  the  plot  of  the 
values  given  for  *°8iokcalc  shows  a  *ar8e  degree  of  scatter.  The  source  of 
this  discrepancy  is  unknown  to  the  present  writer.  Until  it  is  explained  it 
would  probably  be  best  to  disregard  the  data  given  in  Reference  C-3. 

The  group  at  American  Cyanamid^  decomposed  RDX  in  reactors  of 

constant  known  volumes  and  followed  the  reaction  by  gas  chromotographic 
analysis  of  the  undecomposed  RDX.  This  has  an  advantage  over  manometrlc 
determination,  since  it  is  not  subject  to  errors  due  to  secondary  reactions 
among  gaseous  products,  although  the  manipulations  necessary  to  quench  the 
reaction  in  each  reactor  individually  appear  to  have  limited  the  temperature 
range  to  20°.  The  rate  constants  were  independent  of  initial  mass,  reactor 
volume  or  added  pressure  of  NO2  or  helium. 


The  results  described  in  this  series  of  reports  illustrate  the  importanc 
of  statistical  analysis.  The  preliminary  valuesC-8,9  (log  A-19.2,  Ea**48.7) 
were  obtained  from  a  preliminary  analysis,  but  a  subsequent  least  squares 
analysisc-1°  yielded  log  A  -  22.4  and  Ea  •  55.6  kcal/mole.  The  difference 
seems  large  enough  to  be  chemically  significant.  This  illustrates  the 
importance  of  applying  statistical  analysis  when  estimating  Arrhenius 
parameters  from  kinetic  data,  since  it  is  evidently  difficult  to  "eyeball"  a 
slope  and  intercept  with  sufficient  accuracy. 


Joyner''  carried  out  manometrlc  studies  on  the  decomposition  of  RDX; 
the  data  were  believed  satisfactory  despite  some  problems  with  sublimation  and 
pressure  fluctuation.  In  agreement  with  this  the  three  points,  plotted  for 
neat  liquid  RDX  appear  to  fall  on  the  same  line  (main  report.  Figure  1)  and 
the  rate  constants  seem  to  agree  reasonably  well  with  those  of  Robertson1'  ’ 
and  of  the  American  Cyanaroid  grown. C”10 

Thermal  analysis  studies  will  be  considered  in  two  parts  -  first  the 
isothermal  DSC  and  DTA,  then  dynamic  or  temperature-programmed  studies.  We 
begin  with  isothermal  DSC  results.  The  best  isothermal  DSC  result  is  that  of 
Rogers,  who  obtained  log  A  (sec-*)  ■  18.3  and  Efl  -  47.1  kcal/mole  from 
experiments  in  which  correction  was  made  for  gas-phase  reaction  which  proceeds 
concurrently  with  the  liquid  reaction.  To  this  writer's  knowledge,  this  is 
the  only  study  in  which  a  quantitative  correction  for  vapor-phase  reaction  has 
been  parried  out.  Comparison  with  the  uncor.rected  valuesc-*^’*^  of  log  A 

(sec  )  «  16.4  and  E  „  43.1  kcal/mole  reveals  that  the  correction  is  of  the 
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order  of  a  factor  of  100  in  A  and  a  difference  of  4  kcal/mole  in  aC'..  tion 
energy  for  DSC  measurements  of  heat  output  using  the  Perkin-Elme-  ,U0-1B 
apparatus  of  Reference  014;  but  the  magnitude  of  the  correr1' j„  will  of 
course  vary  with  the  mass/volume  ratio  and  the  geometry  of  the  system.  Since 
the  other  isothermal  DSC  studies®"*^*2® .28-30  made  no  attempt  to  correct  their 
results  specifically  for  the  gas-phase  reaction  it  is  not  surprising  to  note 
that  their  activation  energies  and  frequency  factors  fall  somewhat  lower  than 
the  best  DSC®”*^  and  kinetic®"*  ,2,8-10,17  values. 


In  addition  to  the  above  error  (previous  paragraph)  due  to  vapor-phase 
reaction,  non-iso thermal  thermal  analysis  experiments  would  be  affected  by  a 
variation  with  temperature  of  the  proportion  of  the  vapor-phase  reaction;  this 
is  in  addition  to  any  errors  introduced  in  the  derivation  of  the  Kissinger 
equation  (see  p.  17-19)  or  other  method  used  to  estimate  kinetic  parameters. 
None  of  these  studies  attempts  to  correct  for  gas-phase  reaction. 

Consequently,  while  all  non- iso thermal  thermal  analysis  studies  of  which  this 
writer  is  aware  have  been  tabulated  they  will  not  be  discus;  '  individually. 


Activation  energies  for  decomposition  of  RDX  in  solution  (Table  C-2)  seem 
to  be  generally  several  kcal/mole  lower  than  the  values  (Table  Cl)  for  the 
neat  liquid.  Robertson®"*  found  values  for  log  A  (sec"*)  and  activation 
energy  ( ca  15.5  and  41  kcal/mole  respectively)  that  were  noticeably  lower  than 
his  values  (18.5  and  47.5)  for  the  neat  liquid.  However,  the  values  in  TNT 
may  be  affected  by  concomitant  decomposition  of  TNT;  see  following  paragraph. 

The  work  of  Wilby®"22 *2^  was  done  over  a  short  (17-18°)  temperature 
range,  using  large  samples  which  might  have  led  to  self-heating;  furthermore, 
Wilby  pointed  out  that  TNT  was  decomposing  at  the  same  time  as  RDX  and  that 
this  might  have  complicated  matters.  However,  the  Arrhenius  parameters  were 
in  reasonable  agreement  with  those  of  Robertson.®”*  This  agreement  for  such 
large  sample  size  may  suggest  that  self  heating  is  not  a  source  of  error  for 
decomposing  RDX  in  TNT  and  possibly  in  inert  solvents  generally;  however,  note 
the  short  temperature  range  in  Wilby’ s  measurements.  The  Russian  workers  also 
found®"2-*  »2**  that  log  A  (sec"*)  and  Ea  were  lower  in  solution  in 
dinitrobenzene  than  were  the  best  values  for  liquid  RDX. 


The  value  of  '.7.13  obtained  by  Joyner®"*^  is  somewhat  higher  than  the 
value  of  about  14-.6  obtained  in  the  other  entries  of  Table  C-2.  However, 
this  is  to  be  expected  since  it  was  obtained  by  assuming  an  activation  energy 
(44.2  kcal/mole)  determined  for  neat  liquid  RDX;  this  value  is  higher  than 
most  of  those  in  Table  C-2  and  would  thus  be  expected  to  lead  to  a  higher 
frequency  factor. 


The  above  studies  were  all  done  manoraetrically.  Colman  and  Rauch®"2* 
studied  the  decomposition  in  TNT  kinetically,  by  analysis  of  undecomposed  RDX; 
this  would  reduce  errors  due  to  vaporization  followed  by  decomposition  of  RDX 
in  the  vapor  phase.  This  is  especially  pertinent  since  while  RDX 
decomposition  seems  unaffected  by  TNT,  the  TNT  decomposition  is  accelerated  in 
the  presence  of  more  than  a  5%  concentration  of  decomposing  RDX;  the 
composition  of  the  gaseous  products  appears  to  be  a  superposition  of  those 
from  RDX  and  TNT  alone.  Thus,  while  Robertson's  values  for  1%  and  5%  RDX  in 
TNT  are  probably  not  in  error  from  this  source,  studies  using  more  than  about 
5-10%  RDX  and  TNT  may  be  in  error.  However,  the  errors  in  Ea  and  Log  A  may 


not  be  too  great  since  it  was  found^-^  that  activation  energies  for 
decomposition  of  RDX  and  of  solvent  were  about  the  same  within  the  rather 
large  (±  5  kcal/mole)  error  limits. 

Debenham  and  OwenC-^  also  followed  the  reaction  by  analysis  of 
undecomposed  RDX.  Considering  the  short  (10°)  temperature  range,  their 
approximate  activation  energy  of  46  kcal/mole  seems  essentially  the  same  as 
the  other  values  in  Table  C-2  although  a  bit  on  the  high  end  of  the  range. 

Rogers  and  Morris**-11  obtained  an  activation  energy  of  48.4  kcal/mole  by 
DSC,  in  Convalex-10  (polyphenyl  ether).  This  value  is  about  the  same  as  the 
best  values  in  Table  C-2  for  the  decomposition  of  neat  liquid  RDX,  but 
significantly  lower  than  the  value  (67.4  kcal/mole)  reported  on  the  same 
apparatus  by  Rogers  and  Morris**- 11  for  decomposition  of  liquid  RDX;  thus,  the 
activation  energy  for  RDX  in  solution  is  still  less  than  for  the  neat  liquid, 
when  the  two  are  determined  by  the  same  method  on  the  same  apparatus. 

Table  C-3  lists  the  rate  constant  and  temperature  data  used*'-1,1*1,1^,17 
in  what,  on  the  basis  of  the  above  discussion,  the  present  writer  feels  to  be 
the  most  reliable  studies  on  decomposition  of  liquid  RDX.  These  data  were 
analyzed  together  on  a  common  axis  (main  report.  Figure  l,  Table  and 
accompanying  discussion)  to  give  the  best  values  for  log  A  (sec-*)  and  Ea 
reported  there.  Table  C-3  also  gives  the  values  obtained  when  the  programs 
used  in  the  present  report  were  applied  to  the  data  from  the  it.  ividual  data 
sets;  values  reported  by  the  individual  authors  are  given  in  Taole  C-l. 


Table  C-4  lists  the  actual  rate  constant  and  temperature 
data,  used  in  what  the  present  writer  feels,  on  the  basis  of  the 

above  discussion,  to  be  the  most  reliable  data  for  decomposition  of  RDX  in 
solution  in  TNT.  These  data  were  plotted  together  on  a  common  axis  (main 
report.  Figure  5)  and  the  data  were  analyzed  statistically  as  described  in  the 
main  report  (Table  and  accompanying  discussion)  to  give  the  suggested  values 
for  log  A  (sec-1)  and  Fa  reported  there.  Table  C-3  also  gives  the  values 
obtained  when  the  programs  used  in  the  present  report  were  applied  to  the  data 
from  the  individual  data  sets;  values  reported  by  the  individual  authors  are 
given  in  Table  C-2. 
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200  473  5.30  4  Singleton);  ordinary  least  squares  gave  log  A< 

210  483  13. 5  *  47,6  i  1.3  kcal/eole. 

220  493  37.4  * 

230  503  132  * 
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APPF-  “  j. 

REPORTED  ARRHENIUS  PARAMETERS  FOR  THERMAL  DECOMPOSITIONS  OF 
HMX  AND  RDX  AT  TEMPERA^  JRES  BELOW  THEIR  MELTING  POINTS. 

Available  values®”*"^  for  decomposition  of  HMX  and  RDX  at  temperatures 
below  their  melting  points  are  summarized  in  Tables  D-l  and  D-2. 
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As  discussed  elsewhere  (Section  IV.  D.  of  the  main  report),  decomposition 
of  HMX  and  RDX  at  temperatures  below  their  melting  points  is  most  likely  a 
complex  process  involving  the  liquid  and  vapor  states  as  well  as  the  solid 
state.  Furthermore,  the  activation  energies  and  frequency  factors  appear  to 
be  dependent  on  factors  like  degree  of  spreading  of  the  sample  (see  for 
example  Reference  D-7)  which  are  often  difficult  to  evaluate  in  view  of  the 
limited  information  given  in  many  of  the  reports.  Therefore,  although  an 
attempt  has  been  made  to  include  all  available  values  in  Tables  D-l  and  D-2, 
we  will  confine  ourselves  to  pointing  out  some  of  the  main  trends  in  the  data. 

One  trend  is  a  tendency  for  activation  energies  measured  at  or  near  the 
melting  point  by  TGA0-^,  Dgcl>“9 , 10 ,2 1 ,23,2A  or  pj^s^l  techniques  to  have 
what  seem  to  be  unreasonably  high  values,  sometimes  well  over  200  kcal/mole. 
Possibily  these  values  are  artifacts  resulting  from  complexity  in  the  behavior 
of  the  material  as  it  simultaneously  melts  and  decomposes  (see  p.  35  of  the 
main  report). 

Another  pattern  is  that  many®-^ ,14,15,17-20,25  { but  not  all)  of  the 
studies  on  HMX  decomposition  show  changes  in  activation  energy  in  the  region 
around  approximately  245  ±  15°C.  The  magnitude  and  exact  locations  of  these 
changes  are  not  reproducible  from  study  to  study,  except  that  they  generally 
fall  in  the  approximate  range  230-60#C. 

An  attempt0-^®  at  applying  autocatalytic  solid  reaction  kinetics  to  HMX 
decomposition  yielded  an  activation  energy  of  140  kcal/mole  and  a  frequency 
factor  of  4.76  x  10®'’.  However,  the  significance  of  these  numbers  seems 
uncertain;  for  one  thing,  the  k^  in  equation  6.4.2  of  Reference  D-29,  which  is 
equated  by  the  authors  of  Reference  D-28  to  a  rate  constant,  is  described  in 
Reference  D-29  as  an  "adjustable  parameter;"  the  justification  for  relating  it 
to  a  rate  constant  is  not  obvious.  Also,  according  to  this  treatment  the  rate 
constant  comes  out  with  units  of  time  raised  to  the  power  of  -3.5,  although 
rate  constants  generally,  because  of  the  manner  of  their  derivation,  are 
expected  to  be  expressed  in  terms  of  reciprocal  (time  multiplied  by 
( concentration0- 1 ) .  If  the  term  k  tB  in  equation  6.4.2  of  Reference  D-29  is 

replaced  by  (ktl^,  and  k  (-  (k  ) ^  ®  is  considered  as  a  rate  constant  having 
dimensions  sec"  ,  it  would  yield  values  of  Login  A  and  Ea  cliat  would  he  lower 
by  a  factor  of  0  (—3.5,  Reference  D-28)  than  those  (logjg  A  -  35  and  Ea  »  140 

kcal/mole)  given  in  Reference  D-28,  or  logjg  A(sec-*)  »  10  and  Ea  »  40 
kcal/mole. 

Belyaeva  et  ai.0-^  have  reported  a  series  of  manometric  studies  on  HMX 
and  RDX  decomposition  in  which  separate  account  was  taken  of  the  gas  and 


D— 28 

A.  E.  Axworthy,  J.  E.  Flanagan,  and  J.  C.  Gray,  “Interaction  of  Reaction 
Kinetics  and  Nitramine  Combustion,”  Rocketdyne  Division,  Rockwell 
International,  Canoga  Park,  CA,  Report  AFATL-TR-30-58 ,  May  1980  (AD-B052 
86 1L). 

^^M.  Boudart,  “Kinetics  of  Chemical  Processes,”  Prentice-Hall,  New  York, 
Chapter  6,  1968. 
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condensed  phase  decompositions;  the  gas-  and  condensed-phase  rates  were 
evaluated  separately  from  the  following  equation: 

If  -  «c ,1  P2  +  BjKjU  -  n)  ■ 

IB 

O 

where  A  ■  273/760^  T,  B  ■  273R/760M;  K,  and  K2  are  the  decomposition  rates  In 
the  gas  and  condensed  phases  respectively;  V  is  the  volume  of  products  up  to  a 
given  point  in  time;  is  the  final  volume  of  conversion  products;  mo  is  the 
sample  weight,  and  N  is  molecular  weight;  and  dn/dt  is  the  total  rate  of 
thermal  decomposition. 

Initial  rate  data  in  the  gas-  and  condensed  phases  are  given  in  Appendix 
A  and  in  the  final  lines  of  Tables  D-l  and  D-2  respectively.  The  condensed- 
phase  Arrhenius  parameters  in  Tables  D-l  and  D-2  are  based  on  statistical 
analyses  done  by  the  present  writer;  the  original  Russian  authors  of  Reference 
D-22  did  not  report  Arrhenius  parameters  from  their  condensed  phase  data;  it 
should  be  noted  that  this  may  well  be  the  result  of  well-founded  doubts  as  to 
their  exact  significance.  Nevertheless,  it  seemed  worthwhile  to  carry  out  the 
statistical  analyses  whose  results  are  summarized  in  the  last  two  lines  of 
Tables  D-l  and  D-2.  The  “condensed  phase"  Arrhenius  parameters  or  temperature 
coefficients  of  reaction  are  significantly  lower  for  HMX  than  for  RDX;  in  fact 
the  frequency  factor  for  HMX  (log  A(sec"*)  *9.1,  Table  D-2)  is  in  a  range 
that  in  the  gas  phase  would  be  characteristic  of  a  bimolecular  reaction,  while 
that  for  (log  A(sec”*)  ■  13. 5)  is  in  a  region  that  in  the  gas  phase  would  be 
more  characteristic  of  a  unimolecular  reaction.  Possible  reasons  for  this 
discrepancy  Include  (a)  differences  In  chemical  mechanisms  for  the  solid-state 
decomposition  of  HMX  and  RDX;  and  (b)  differences  in  kind  or  contribution  of 
physical  mechanisms,  in  particular  vaporization  and  gas-phase  decomposition  of 
KMX  and  RDX,  and  liquefaction  of  the  solid  phase  due  to  condensation  of 
products  such  as  hydroxymethylformamlde. 


3 


f- 

< 

X 


§ 


b. 

C 


-  x~s  O 


.1 

o  c 

^  5 

M  o 


Al? 


o 

c 


TJ 

CO 

t 

u 

a 

<0 

c 

3 


6C 


6C 

CM 

C 

<N 

z 


•o 

CO 

c 

ix 

a 

w 


CN 

C 


■o 

CO 

a 

u 

a. 

io 

c 

3 


U 

CN 


C 


*3 

CO 

a 

u 

a. 

«5 

C 


41  41 
Cl  Cl 
£  £ 
«0  CO 
CO  to 


6C 


CN  C*> 

A  A 


ii 


aj 

:r 


x 

i* 

y 

4i 

E  V  V 

u  C  «o  <0 

C  u  -i  -« 

-  *-  U  lx 

u  u 

CO  at  41  D 

V  Cl  Id  U 

W*  cfl  3  3 

(0  VJ 

3  T.  V>  * 

O  cn  ai  41 

xX  nj  Id  L. 

u.  6  a.  £ 


A 

*cr 

A 


4) 

U 

3 


cn 

41 

u 

a. 


A 


A 


r-> 

A 


A 


O 

CN 


:  Cl  c 
>.  -* 
H  X 


41 

o 

CJ 

< 


a*  ii 

o.  *J  ^ 
SCO 


102 


180-98  453-71  Maximum  (p  max)  63.6  266  Pressure  rise  D-7  Sample  spread,  100  mm  N2»  0.2 

sample 


130-80  403-33  InltlaKup  to  O.IX)  11.?  38.3  165  Proaouro  rise  D-22 


'A 


G 

& 


8 

u 

m 


$ 

X 

a 


a d 

E 

X 

£ 


W 

J 

3 


•4  « 

J! 
S  S 

■“i 


p-  ■ 

§  K 


©  3 
«*-  o 

s  °r 

Sis 


I  * 
-  H 
12  J 
!£u‘ 


? 
J 

>  A  o 


0 

« 

0  <J 

-4  C 
<0  V 
Sn  U 
^  © 

22 


«  S 

13 


IM  |  >© 

TJ  — • 
<U  -4 
£  — * 

4-*  O  e~ 

•  X 

£  u  2 

3 

V  U  •*- 

»  o 

0 

J=  «M  C 

«  &  e 

■*4 

•«  m  w 
e  «i  **4 

<u  .  j  « 

•  o 

®  a  a 

«  <H  b 

3  4-1  O 


V  *>» 

62 


>  T>  • 

©  c  «  Jp 

«  «  4J 

©  o  «j  o 

£  «  «  S 

U  -O  ■  N 


3 


X 

u  e 

•  o 


$2 


|4  «  C 
M  U  O 

I  2S 

O  u  O 

Jo  m 

2£ 


104 


n 


t Ion  only;  50  a*  simple 


245-70  518-543 


222”*2  *95-515  42  175  Haas  Spec.  Study 

of  decomposition  products 

251-64  524-37  42*2  176+8  Isothermal  DSC  D-27 


/  • 


APPENDIX  E 


COMPUTER  PROGRAM  USED  FOR  LEAST-SQUARES  FIT  OF  ORIGINAL  DATA 

PROGRAM  ACTEN 

V)  I  MENS  I DI4  TITL  E  (BO) 

DH  30  III-- 1,10 

HE  AD ( 5 , 200 1 , END=5000 , EER-5000  >  T I TLE 
, •«»<•>  1  FORMAT  (SOAl  ) 

WE  1 TE ( 6 , 2002 ) T I T  LE 
200 2  FORMAT  (BOA  I  ,  /  ) 

I J--0 
RX-O. 

GX2--0. 

SV  =  0 . 

Sr  2--0. 

3  X  Y  U  - 

I'd  10  I  ~  1,1 000 

READ  ( 5 ,  ►  ,  END*=5000 ,  ERR=»5000 )  X  ,  V 
•«. "  i')  FORMAT  < 2F  10.4) 

IT  (X.EQ.999. )G0  TO  20 
WR I TE 16,2000) X , Y 
N-N+  1 
SX-SX+X 
GX2«SX2+X**2 
.  9V~bV ■+ V 
BY2--BY2+Y*  *2 
GXY-SXY+X*Y 
h-  CONTINUE 
20  CONTINUE 

I  1-SX2~< (SX**2) /FLOAT (N) ) 

I  '2- CSXY- C  < SX  *SY > /FLOAT <N> ) ) /FI 
F 3-  <  -~S  X *F2+SY )  /FLOAT  <N) 

M=<  <F?**2*SX2)  -  <SXY*2*F2)  +SY2-FL0AT  <N)*F3**2>  /FLOAT  <N> 

FS---F4/F  1 

EH -SORT (SX2*F5/FL0AT (N) ) 

D--F3 

S/M=BQET <F5> 

«  IT— F  2 

WE  IT  E 1 6 , 400 ) 

WRITE (6,200) XM , SXM , B , SO 
AC  T  F  AC-2 . 3026# 1 . 9072 
A/M---  ( - 1 .  )  *XM*ACTFAC. 

ASXfv=S  XM *  ACT  FAC 
WE  t  TE ( 6 , 502 ) AXM , ASXM 
30  CONTINUE 
200  I  OEM A f ( 4F 1 1.4) 

300  FORMAT (/) 

-I  «.>U  FORMAT  </,'  SLOPE,  SIS  SLOPE,  INTERCEPT,  SIG  INTERCEPT 
302  FORMA T</,'  ACT  EN  SIG  ACT  EN ' , / / , 2F 10. 4 , / / ) 

50 uO  CONTINUE 
END 
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